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Abstract 

We establish sharp energy decay rates for a large class of nonlinearly first-order damped 
systems, and we design discretization schemes that inherit of the same energy decay rates, uni¬ 
formly with respect to the space and/or time discretization parameters, by adding appropriate 
numerical viscosity terms. Our main arguments use the optimal-weight convexity method and 
uniform observability inequalities with respect to the discretization parameters. We establish 
our results, first in the continuous setting, then for space semi-discrete models, and then for 
time semi-discrete models. The full discretization is inferred from the previous results. 

Our results cover, for instance, the Schrodinger equation with nonlinear damping, the 
nonlinear wave equation, the nonlinear plate equation, the nonlinear transport equation, as 
well as certain classes of equations with nonlocal terms. 

Keywords: stabilization, dissipative systems, space/time discretization, optimal weight convexity 
method. 
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1 Introduction 

Let X be a Hilbert space. Throughout the paper, we denote by || • ||x the norm on X and by (•, • )x 
the corresponding scalar product. Let A : D(A) ->lbea densely defined skew-adjoint operator, 
and let B : X —> X be a nontrivial bounded selfadjoint nonnegative operator. Let F : X —> X be 
a (nonlinear) mapping, assumed to be Lipschitz continuous on bounded subsets of X. We consider 
the differential system 

u'(t) + Au(t) + BF(u(t)) = 0. (1) 

If F = 0 then the system (1) is conservative, and for every uq £ D(A), there exists a unique 
solution u(-) € C°(0, +oo; D(A)) fl C 1 (0, +oo; X) such that u(0) = uq, which satisfies moreover 
IKOIU = IK0 )IIa-, for every t > 0. 

If F ^ 0 then the system (1) is expected to be dissipative if the nonlinearity F has “the good 
sign". Defining the energy of a solution u of (1) by 

Eu(t) = l\Ht)\\x, ( 2 ) 

we have, as long as the solution is well defined, 

E' u (t)=-(u(t),BF(um x . (3) 

In the sequel, we will make appropriate assumptions on B and on F ensuring that E' u (t) ^ 
0. It is then expected that the solutions are globally well defined and that their energy decays 
asymptotically to 0 as t —> +oo. 

The objective of this paper is twofold. 

First of all, in Section 2 we provide adequate assumptions under which the solutions of (1) have 
their norm decaying asymptotically to 0 in a quasi-optimal way. This first result, settled in an 
abstract continuous setting, extends former results of [9] (established for damped wave equations) to 
more general equations and damping operators for stabilization issues based on indirect arguments 
(for direct arguments see e.g. [5, 8]). 

Then, in Section 3, we investigate discretization issues, with the objective of proving that, for 
appropriate discretization schemes, the discrete approximate solutions have a uniform decay. In 
Section 3.1, we first consider spatial semi-discrete approximation schemes, and in Section 3.2 we 
deal with time semi-discretizations. The full discretization is done in Section 3.3. In all cases, we 
establish uniform asymptotic decay with respect to the mesh size, by adding adequate viscosity 
terms in the approximation schemes. 
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2 Continuous setting 

2.1 Main result 

Assumptions and notations. First of all, we assume that 


{u, BF(u))x ^ 0, 


( 4 ) 


for every u € X. Using (3), this first assumption ensures that the energy E u (t ) defined by (2) is 
nonincreasing. 

Since B is bounded, nonnegative and selfadjoint on X , it follows from the well-known spectral 
theorem that B is unitarily equivalent to a multiplication (see, e.g., [23]). More precisely, there 
exist a probability space Q with measure p, a real-valued bounded nonnegative measurable function 
b defined on f2 (satisfying ||&||,L~(n,/q = ||-B||), and an isometry U from L 2 (VL,p) into X , such that 

(U~ 1 BU f)(x) = b(x)f(x), 


for almost every x £ f2 and for every / £ p). Another usual way of writing B is 

r+oo 

/ 0 


B = 


XdE(X), 


where the family of E( A) is the family of spectral projections associated with B. We recall that 
the spectral projections are obtained as follows. Defining the orthogonal projection operator Q\ 
on L 2 (n,p) by ( Q\f)(x ) = for every / £ L 2 (fl,p) and for every i £ SI, we have 

E( X) = UQxU- 1 . 

We define the (nonlinear) mapping p : L 2 (Q,p) —>- L 2 (VL,p) by 

p(f) = U~ 1 F(U /), 


for every / £ L 2 (Q,p). In other words, the mapping p is equal to the mapping F viewed through 
the isometry U. 

Note that, setting / = U~ 1 u, the equation (1) is equivalent to dtf + Af + bp(f) = 0, with 
A = U~ 1 AU, densely defined skew-adjoint operator on U 2 (fl,/z), of domain U~ 1 D{A). 

We assume that p(0) = 0 and that 

fp(f) > 0, (5) 

for every / £ L 2 (U, p). Following [5, 8, 9], we assume that there exist ci > 0 and C 2 > 0 such that, 
for every / £ L°°{Q,p), 

ci 3(1/(201) < \p(f)(%)\ < c 2 g _1 (|/(x)|) for almost every x £ U such that |/(x)| < 1, 

(o) 

Ci |/(x)| ^ |jo(/)(ar)| < c 2 |/(a;)| for almost every x £ U such that |/(x)| ^ 1, 


where g is an increasing odd function of class C 1 such that 

9 ( 0 ) = «'( 0 )= 0 , 

ff(s) 

and such that the function H defined by F[(s) = y/sg{y/s), 
on [0, Sq] for some sq £ (0,1] (chosen such that g(so) < 1). 


- > 0 

s —^0 

for every s £ [0,1], is strictly convex 


3 




We define the function H on ]R by H(s) = H(s) for every s G [0,Sg] and by H(s) = +oo 


otherwise. Using the convex conjugate function H* of H , we define the function L on [0, +oo) by 
L{ 0) = 0 and, for r > 0, by 



(7) 


By construction, the function L : [0, +oo) —> [0, Sg) is continuous and increasing. We define the 
function A h ■ (0, s§] —> (0,+oo) by A h(s) = H(s)/sH'(s), and for s ^ 1 /H'(sq) we set 



( 8 ) 


The function if : [\/H'(s(f), +oo) — > [0,+oo) is continuous and increasing. 

Throughout the paper, we use the notations < and ~ in many estimates, with the following 


meaning. Let S be a set, and let F and G be nonnegative functions defined on 1R x fi x 5. The 


notation F < G (equivalently, G > F) means that there exists a constant C > 0, only depending on 
the function g or on the mapping p, such that F(t, x, A) ^ CG[t , x, A) for all (t,x, A) SlRxfix S. 
The notation F\ ~ F -2 means that F\ < F 2 and F\ > F?. 

In the sequel, we choose S = X, or equivalently, using the isometry U, we choose S = L 2 (U, p), 
so that the notation < designates an estimate in which the constant does not depend on u £ X, 
or on / e L 2 (U, p), but depends only on the mapping p. We will use these notations to provide 
estimates on the solutions it(-) of ( 1 ), meaning that the constants in the estimates do not depend 
on the solutions. 

For instance, the inequalities ( 6 ) can be written as 


5(1/1) 1$ \p(f)\ 9 1 (|/|) on the set I/I < 1, 

Ip(/)I — I/I on the set I/I > 1 . 

The main result of this section is the following. 


Theorem 1 . In addition to the above assumptions, we assume that there exist T > 0 and Ct > 0 
such that 



for every solution of <f>' (t)+Acf>(t) = 0 (observability inequality for the linear conservative equation). 

Then, for every u$ G X, there exists a unique solution u (-) G C°( 0, + 00 ; J A)nC 1 (0, + 00 ; D(A) r ) 
of ( 1 ) such that u(0) = mq . 1 Moreover, the energy of any solution satisfies 



( 10 ) 


for every time t ^ 0, with 71 ~ ||B ||/72 and 72 ~ Ct/(T' 3 ||.B 1 / 2 || 4 + T). If moreover 


limsupAij(s) < 1, 


( 11 ) 


then we have the simplified decay rate 


E u (t) < T max( 7 l ,U u (0)) (iT )" 1 , 


for every time t > 0 , for some positive constant 73 ~ 1 . 

1 Here, the solution is understood in the weak sense, see [13, 17], and D(A)' is the dual of D(A) with respect to 

the pivot space X. If uq G D(A ), then u(-) E C°(0, +oo;Z)(A)) D C 1 (0, + 00 ; X). 
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Theorem 1 improves and generalizes to a wide class of equations the main result of [9] in 
which the authors dealt with locally damped wave equations. The case of boundary dampings 
is also treated in [8] (see also [5]) by a direct method, which provides the same energy decay 
rates. The result gives the sharp and general decay rate L(l/'0 _ 1 (t)) of the energy at infinity 
(forgetting about the constants), and the simplified decay rate (H')~ 1 ( 1/t) under the condition 
(11). It is also proved in [ 8 ] that, under this condition, the resulting decay rate is optimal in 
the finite dimensional case, and for semi-discretized nonlinear wave or plate equations. Hence 
our estimates are sharp and they are expected to be optimal in the infinite dimensional case. 
The proof of optimality relies on the derivation of a one-step decay formula for general damping 
nonlinearity, on a lower estimate based on an energy comparison principle, and on a comparison 
lemma between time pointwise estimates (such our upper estimate) and lower estimates which are 
of energy type. Note also that p has a linear growth when g'(0) ^ 0. In this case the energy 
decays exponentially at infinity. Moreover, even in the finite dimensional case, optimality cannot 
be expected when limsupA#(s) = 1 (functions p leading to that condition are close to a linear 

growth in a neighborhood of 0 ), and it is possible to design examples (linear feedback case) with 
two branches of solutions that decay exponentially at infinity but do not have the same asymptotic 
behavior. 

Let us recall some previous well-known results of the literature. The first examples of nonlinear 
feedbacks were only concerning feedback functions having a polynomial growth in a neighborhood 
of 0 (see e.g. [36, 24] and the references therein). As far as we know, the first paper considering 
the case of arbitrary growing feedbacks (in a neighborhood of 0) is [26]. In this paper, the analysis 
is based on the existence (always true) of a concave function h satisfying h(sp(s )) ^ s 2 + p 2 (s) for 
all |s| ^ N (see (1.3) in [26]). The paper is very interesting but provides only two examples of 
construction of such function h in Corollary 2, namely the linear and polynomial growing feedbacks. 
The results use only the Jensen’s inequality (not the Young’s inequality), and allow the authors to 
compare the decay of the energy with the decay of the solution of an ordinary differential equation 
S'(t) + q(S(t)) = 0 where q(x) = x — (/ + p)~ l (x) and p(x) = (cl + h(Cx))~ 1 (Kx) where c, C , K 
are non explicit constants and / -1 stands for the inverse function of /. In the general case, these 
results do not give the ways to build an explicite concave function satisfying h(sp(s)) ^ s 2 + p 2 (s). 
No general energy decay rates are given in an explicit, simple and general formula, which besides 
this, could be shown to be "optimal". Due to this lack of explicit examples of decay rates for 
arbitrary growing feedbacks in other situations than the linear or polynomial cases, other results 
were obtained, also based on convexity arguments but through other constructions in [32, 33] 
(see also [37]) through linear energy integral inequalities and in [30] through the comparison with 
a dissipative ordinary differential inequality. In both cases, optimality is not guaranteed. In 
particular, [32, 33] do not allow to recover the well-known expected "optimal" energy decay rates 
in the case of polynomially growing feedbacks. Optimality can be shown in particular geometrical 
situations, in one dimension when the feedback is very weak (as for p(s) = e _1 ' s for s > 0 close to 
0 for instance), see e.g. [42, 5]. Hence the challenging questions are not only to derive energy decay 
rates for arbitrary growing feedbacks, but to determine whether if these decay rates are optimal, 
at least in finite dimensions and in some situations in the infinite dimensional case, and also to 
derive one-step, simple and semi-explicite formula which are valid in the general case. This is the 
main contribution of [5, 8 ] for direct methods and of the present paper for indirect methods for 
the continuous as well as the discretized settings (see also [9] for the continuous setting). Note also 
that the direct method is valid for bounded as well as unbounded feedback operators. 

Several examples of functions g (see ( 6 )) are given in Table 1, with the corresponding Ah and 
decay rates. Note that other examples can be easily built, since the optimal-weight convexity 
method gives a general and somehow simple way to derive quasi-optimal energy decay rates. 
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9(s ) 

A H (s) 

decay of E(t) 

s/ ln -p (l/s), p > 0 

limsupA#(s) = 1 

e -t 1/{p+1> / t i/( P +U 

s p on [0, Sq]> P > 1 

A h(s) = < 1 

t -2/(p-l) 

e- 1 ^ 

lim A h(s) = 0 

s \,0 

V Ht) 

s p ln 9 (l/s), p > 1 , q > 0 


t -2/(p-l) ln -2<J/(P-I)( t ) 

e~ in p (l/s) , p > 2 

lim A h(s) = 0 

s \0 

e - 21 n 1 /P(i) 


Table 1: Examples 


In the four last examples listed in Table 1, the resulting decay rates are optimal in finite 
dimension and for the semi-discretized wave and plate equations. Moreover, (11) is satisfied. 

Remark 1. It also is natural to wonder whether the linear feedback case, that is g(s) = s for every 
s close to 0 , as well as the nonlinear feedback cases which have a linear growth close to 0 , such as 
the example g(s) = arctans for s close to 0 , are covered by our approach in a "natural continuous" 
way. These cases can be treated under a common assumption, which is indeed g'( 0) 0. Notice 

that an apparent difficulty lies in the fact that the function H is the identity function and is not 
strictly convex anymore, so that our construction may seem to fail. This limit case is investigated 
in the following result, whose proof is postponed at the end of Section 2.3. It is obtained under 
slight modifications in the proof of Theorem 1. This approach is also valid for the direct approach 
presented in [5, 8], which leads to continuous nonlinear integral inequalities. We will also formulate 
a general result in this direction below. 

Corollary 1 . Let us assume that g'( 0) ^ 0. Under the same assumptions as those of Theorem 1 
(namely, the ones at the beginning of Section 2.1 as well as the observability inequality (9 )), there 
holds limsups^Q Ah (s) = 1 and 


E u (t) < T max( 7 i,£' tt ( 0 )) exp(~ 7 2 t), 
for every time t > 0 with 71 ~ ||R ||/72 and 72 ~ Ct , /(T 3 ||I ? 1 / 2 || 4 + T). 

Let us now apply this generalization to the direct optimal-weight convexity method as intro¬ 
duced in [5, 8 ]. In all the results presented in these two papers, and when the bounded as well as 
unbounded feedback operator p satisfies ( 6 ) with a function g such that g( 0 ) = 0 = g'(0) we can 
extend the given proofs to the cases for which g{ 0 ) = 0 whereas g'(0) 0 (i.e. when g has a linear 

growth close to 0). More precisely, we can extend the proof of Theorems 4.8, 4.9, 4.10 and 4.11 
in [8] (but also in the more general framework as presented in Theorem 4.1) to the case for which 
g\ 0) ^ 0. For this, it is sufficient as for the proof of Corollary 1 to replace g by the sequence 
of functions g e defined by g £ (s) = s 1+e for every |s| ^ 1, where e £ (0,1). One can then apply 
the optimal-weight convexity method to this sequence, and define the associate optimal-weight 
function w E (-) = L~ 1 { — ) in a suitable interval (see the above references for more details). We 


then prove that 


2/3 


w e (E(t))E(t)dt ^ ME(S) V 0 < 3 ^ T, 


where M, /3 can be chosen independently on e. We then let e goes to 0. Thanks to the proof of 
Corollary 1, we know that the sequence w e converges poinwise on (0,/3r(j) towards 1. This leads 
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to the inequality 

[ E{t)dt < ME(S) VO < 5 < T, 

Js 

from which we deduce that E decays exponentially at infinity (see e.g. [24]). 

We next provide some typical examples of situations covered by Theorem 1. 

2.2 Examples 

2.2.1 Schrodinger equation with nonlinear damping 

Our first typical example is the Schrodinger equation with nonlinear damping (nonlinear absorp¬ 
tion) 

idtu(t , x) + A u(t, x ) T ib(x)u(t, x)p{x, \u(t, x)|) =0, 

in a Lipschitz bounded subset O of IR n , with Dirichlet boundary conditions. In that case, we 
have X = L 2 (0, C), and the operator A = —iA is the Schrodinger operator defined on D(A) = 
l?g(0,C). The operator B is defined by ( Bu)(x ) = b[x)u{x) where b £ is a nontrivial 

nonnegative function, and the mapping F is defined by (F(u))(x) = u(x)p(x,\u(x)\), where p is 
a real-valued continuous function defined on Cl x [ 0 ,+oo) such that p(-, 0 ) = 0 on O, p(x,s ) ^ 0 
on Cl x [0, Too), and such that there exist a function g £ C 1 ([—1,1], It) satisfying all assumptions 
listed in Section 2.1, and constants ci > 0 and C 2 > 0 such that 

cig(s) < sp(x, s ) ^ c 2 3 _1 (s) if 0 < s < 1 , 
cis ^ p(x, s ) ^ C 2 S if s ^ 1 , 

for every x £ Q. Here, the energy of a solution u is given by E u (t) = i \u(t, x)\ 2 dx, and we have 
E' u (t) = — fn b(x)\u(t, x)\ 2 p{x, | u(t, x)|) dx < 0. Note that, in nonlinear optics, the energy E u (t) is 
called the power of u. 

As concerns the observability assumption (9), it is well known that, if6(-)^a>0on some 
open subset w of H, and if there exists T such that the pair (co,T) satisfies the Geometric Control 
Condition, then there exists Ct > 0 such that 

Cr||#V)ll! W) J n b(x)\4>(t, x )| 2 dxdt , 

for every solution <f> of the linear conservative equation dt4> — iA(f> = 0 with Dirichlet boundary 
conditions (see [28]). 

2.2.2 Wave equation with nonlinear damping 

We consider the wave equation with nonlinear damping 

dttu(t, x) — A u(t, x) T b(x)p{x, dtu(t , x)) = 0, 

in a C 2 bounded subset f 1 of 1R", with Dirichlet boundary conditions. This equation can be written 
as a first-order equation of the form (1), with X = Hq(D) x L 2 (H) and 



defined on D(A) = Uq(H) fl x Hg(Sl). The operator B is defined by (B(u,v))(x) = 

(0, b(x)v(x)) T where b £ L°°(fl) is a nontrivial nonnegative function, and the mapping F is defined 
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by (F(u,v))(x) = ( 0 , p(x, c(:r))) T , where p is a real-valued continuous function defined on Cl x 1 R 
such that p(-, 0 ) = 0 on Cl, sp(x,s) ^ 0 on Cl x R, and such that there exist a function g £ 
C 1 ([— 1,1], H) satisfying all assumptions listed in Section 2.1, and constants ci > 0 and c 2 > 0 
such that 


c iS(l s l) ^ \p(. x > s )\ ^ c 2ff 1 (l s l) ^ 

ci|s| < | p(x, s)| ^ c 2 |s| if s > 1 , 

for every a; £ fl. Here, the energy of a solution u is given by E u (t) = i f n ((|Vu(t, x )| 2 + (d t u(t, x)) 2 ) dx, 
and we have E' u (t) = — f n b(x)d t u(t, x)p(x, d t u(t, x)) dx ^ 0. 

The framework of this example is the one of [9]. 

As concerns the observability assumption (9), it is well known that, if b(-) ^ a > 0 on some 
open subset u> of fl, and if there exists T such that the pair (co,T) satisfies the Geometric Control 
Condition, then there exists Ct > 0 such that 

Ct\\{4>(0, •)> •))llffi ( n )X L 2 ( fi ) J^b{x)(d t <j)(t,x)) 2 dxdt, 

for every solution </> of the linear conservative equation dtt4> — A <f> = 0 with Dirichlet boundary 
conditions (see [ 12 ]). 

2.2.3 Plate equation with nonlinear damping 

We consider the nonlinear plate equation 

d t tu(t, x) + A 2 u(t, x) + b{x)p{x, d t u(t, x)) = 0, 

in a C 4 bounded subset f l of 11", with Dirichlet and Neumann boundary conditions. This equation 
can be written as a first-order equation of the form (1), with X = ifg(fl) x L 2 (fl) and 



defined on D(A) = (if 2 (f2) fl H 4 (fl)') x Ug(fl). The operator B is defined by (B(u,v))(x) = 
(0, b(x)v(x)) T where b £ L°°( fl) is a nontrivial nonnegative function, and the mapping F is 
defined by (F(u,v))(x) = (0, p(x, u(x))) T , where p is a real-valued continuous function defined 
on Cl x [0, Too) such that p(-,0) = 0 on fl, sp(x,s) ^ 0 on Cl x It, and such that there ex¬ 
ist a function g £ C 1 ([— 1,1], 1R) satisfying all assumptions listed in Section 2.1, and constants 
Ci > 0 and c 2 > 0 such that (12) holds. Here, the energy of a solution u is given by E u (t) = 
i f n ((A u(t,x)) 2 + (d t u(t,x)) 2 ) dx, and we have E' u {t) = — f Q b(x)d t u(t, x)p(x, d t u(t, x)) dx ^ 0. 

The framework of this example is the one of [ 6 ]. 

A sufficient condition obtained in [28], ensuring the observability assumption (9), is the follow¬ 
ing: if &(•) ^ a > 0 on some open subset u of O for which there exists T such that the pair (w,T) 
satisfies the Geometric Control Condition, then there exists Ct > 0 such that 

Ct||(0(O, •), d t c/)(0, •))llff 2 (n)xi2(n) < f f b{x){d t (l>{t,x)) 2 dxdt, 

for every solution cj> of the linear conservative equation duty + A 2 </> = 0 associated to the corre¬ 
sponding boundary conditions. 


2.2.4 Transport equation with nonlinear damping 

We consider the one-dimensional transport equation 

dtu(t, x ) + d x u{t , x) + b(x)p(x, u(t, x)) = 0 , x £ ( 0 , 1 ), 

with periodicity conditions u(t, 0) = u(t, 1). This equation can be written as a first-order equation 
of the form (1), with X = L 2 { 0,1) and A = d x defined on D(A) = {u £ f? 1 (0,1) | u( 0) = u(l)}. 
We make on p the same assumptions as before. The energy is given by E u (t) = \ fg u(t,x) 2 dx, 
and we have E' u (t) = — b(x)u(t, x)p(x, u(t, x)) dx ^ 0. The observability inequality for the 
conservative equation is satisfied as soon as the observability time is chosen large enough. 

On this example, we note two things. 

First of all, the above example can be easily extended in multi-D on the torus T” = M n /Z™, by 
considering the following non-linear transport equation 

d t u(t, x) + di v(v(x)u(t, x)) + b{x)p{x, u(t , x)) = 0, t > 0, x £ T”, 

where v is a regular vector field on T” such that div(u) = 0. The divergence-free condition on the 
function v ensures that the operator A defined by Az = div(v(:r)z(:r)) on 

D(A) = {z£H\ r) I *(■ +e0 =*(■), Vie [l,n]}, 

where ej denotes the i -th vector of the canonical basis of M™, is skew-adjoint. 

Second, in ID we can drop the assumption of zero divergence, by using a simple change of 
variable, which goes as follows. We consider the equation 

d t u(t, x) + v(x)d x u(t, x) + bp(x, u(t, x)) =0, t > 0, x £ (0, 1 ), 

with v a measurable function on (0,1) such that 0 < V- ^ v(x) ^ v+. Then, using the change of 
variable x i —> fg we immediately reduce this equation to the case where v = 1 . 

Hence our results can as well be applied to those cases. 

2.2.5 Dissipative equations with nonlocal terms 

In the three previous examples, the term b(-)p(-,-) is a viscous damping which is local. In other 
words, the value at x of the function F(u) does only depend on the value at x of the function u. 
To illustrate the the potential of our approach and the large family of nonlinearities that it covers, 
We slightly modify here the examples presented in the sections 2.2.1, 2.2.2, 2.2.3 and 2.2.4, by 
providing several examples of viscous damping terms containing a non-local term. 

We refer to the previous sections for the precisions on the boundary conditions and the func¬ 
tional setting associated to each system. 

We consider the non-linear systems 

idtu(t, x) + A u(t, x) + ib{x)u{t, x)p(\u\)(t, x) = 0 
d t tu(t, x ) — A u(t, x) + b(x)p(d t u)(t, x) = 0 
d t tu(t, x) + A 2 u(t, x) + b(x)p(dtu)(t, x) = 0 
dtu(t, x) + d x u(t , x) + b(x)p(u)(t, x) = 0 , 

where the non-linear term p is defined by 

P{f)(x) = <p(f(x),Af(f)), 
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where p : IR 2 —> IR is a continuous function and A f : L 2 (Q) —> IR stands for a non-local term. 
We can typically choose A/"(/) = J n x( x ).f( x ) dx with x € L 2 (Q) whenever fl is bounded or \ 
smooth with compact support else. In the framework of Section 2.2.4, one is also allowed to choose 
A/"(/) = K * / with K £ L 2 { T"). We also impose that f satisfies the following uniform Lipschitz 
property: there exists C > 0 such that 

|<p(s, r) - f(s', r)| + r) - p(s, r')| < C(|s - s'| + |s|.|r - r'|) 

for every (s,s '£ IR 4 . As a consequence, one easily infers that the mapping p is Lipschitz 
from L 2 (fl) into L 2 (fl). 

Moreover, we choose the function p odd with respect to its first variable and such that p(s : r) ^ 
0 for every s ^ 0 and r £ IR. It follows that the assumption fp(f ) ^ 0 is satisfied by every 
/ £ L 2 (fl). 

Finally, we assume that the assumption (6) is satisfied. 

Let us provide an example of such a function p. Notice that, if there exist two positive constants 
k\ and &2 such that for every r £ IR, there exist two positive real numbers c T and C T in [fci,^] 
such that 

p(s , r) ~ c T s 3 as s —» 0 and p(s , r) ~ C T s as s —> +oo, 

then the assumption (6) is satisfied. A possible function p is given by 

p(s, t) = p±(s)p 2 (t) where p\ : IR 3 s >-> s — sin s 

and p 2 : IR —t IR denotes any function bounded above and below by some positive constants, for 
instance p 2 (r) = tt + arctan(r). 

2.3 Proof of Theorem 1 

First of all, note that the global well-posedness follows from usual a priori arguments. Indeed, in 
in sequel we are going to consider the solution, as long as it is well defined, and establish energy 
estimates. Since we prove that the energy (which is the Hilbert norm of A) is decreasing, the global 
existence of weak and then strong solutions follows (see, e.g., [13, Theorem 4.3.4 and Proposition 
4.3.9]). Hence, in the sequel, without taking care, we do as if the solution were globally well defined. 
Note that uniqueness follows from the assumption that F is locally Lipschitz on bounded sets. 
The proof goes in four steps. 

First step. Comparison of the nonlinear equation with the linear damped model. 

In this first step, we are going to compare the nonlinear equation (1) with its linear damped 
counterpart 

z'(t) + Az(t) + Bz(t ) = 0. (13) 

Lemma 1. For every solution u(-) of (1), the solution of (13) such that z( 0) = u(0) satisfies 

WB^zmx dt ^ 2 [ (\\B 1 / 2 u(t)\\x + ||I? 1/2 FHf))|| 2 ,) dt. (14) 

Proof. Setting = u(t ) — z(t), we have 

(ip'(t) + Ai/j(t) + BF(u(t )) - Bz{t), if{t))x = 0. 

Denoting E^(t) = |||V , (^)llx) it follows that 

E'^(t) + ||£ 1 / 2 z(f)|| 2 Y = -(u(t),BF(u(t)))x + {B 1 / 2 F(u(t)), B 1 / 2 z(t)) x + (.B 1 ' 2 u{t),B 1 / 2 z{f)) x . 
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Using (4), we have (uft), BF(u(t)))x ^ 0 , and hence 

E'^i) + \\B^ 2 z(t)\\x < \\B^F{um\x\\B^z(t)\\ x + \\B^ 2 u{t)\\ x \\B^ 2 z(t)\\ x . 

2 ,2 

Using the Young inequality ab < with 6 = ±, we get 

U;(t) + ||S 1 ^ 2 -(t)||jc < ^||-B 1/2 z(t)|| 2 Y + ||S 1 / 2 F(u(t))|| 2 f + ||B 1 / 2 u(t)|| 2 Y , 

and thus, 

E'^t) + 111^/2.^11^ ^ \\B l / 2 F(u(t))\\x + ||B 1/2 u(i)|&. 

Integrating in time, and noting that E^(0) = 0, we infer that 

E AT) + \f o ||i? 1/2 *(t)llx dt < jT (llB^Fiumix + \\B 1/2 u(t)\\x) dt. 

Since E^{T) ^ 0, the conclusion follows. □ 


Second step. Comparison of the linear damped equation with the conservative linear equation. 
We now consider the conservative linear equation 

Aft) + Acfft) = 0 . (15) 

Lemma 2. For every solution z(-) of (13), the solution of (15) such that cf>{ 0) = z(0) is such that 

[ \\B 1/2 At)\\ 2 x dt < k T [ \\B 1/2 z(t)\\ 2 x dt, (16) 

Jo Jo 

with k T = 8 T 2 ||S 1 / 2 || 4 + 2 . 

Proof. Setting Oft) = <f>ff) — z(t), we have 

{O'{t) + A0 ft) - Bz(t), Oft)) x = 0. 

Denoting Eg ft) = 4||0(i)ll\'> it follows that E' g ft) = (Bzft),0ff)) x ■ Integrating a first time over 
[0, t], and a second time over [0, T], and noting that Eg(0) = 0, we get 

[ Egft) dt — f f (Bz(s),0(s))x dsdt = f (T - t)(Bzft),0ft)) x dt. 

Jo Jo Jo Jo 

Applying as in the proof of Lemma 1 the Young inequality with 0 = 1 yields 


Wmfxdt^ T 2 \\Bz(t)\\ 2 x dt + 


1 


wmwxdt, 


10 


and therefore, since B is bounded, 

r T 


\j o WOfOfxdt^WB^f J WB^zmxdt. 

Now, since (fft) = Off) + zft), it follows that 

rj~i rj~t rj~t rj-\ 

[ \\B 1 / 2 cf,(t)\\ 2 x dt^2 [ \\B 1 / 2 0(t)\\xdt + 2 [ \\B 1 t 2 z{f)\\ 2 x dt < k T f \\B^ 2 z{f)\\ 2 x dt. 

Jo Jo Jo Jo 

The lemma is proved. □ 
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Third step. Nonlinear energy estimate. 

Let j3 > 0 (to be chosen large enough, later). Following the optimal weight convexity method 
of [5, 8], we define the function 

= (17) 

for every s G [0,/?Sq). In the sequel, the function w is a weight in the estimates, instrumental in 
order to derive our result. 

Lemma 3. For every solution u(-) of (1), we have 

£w(E^m (\\B^u{t)\\ 2 x + WB^Fiumix) dt 

<T\\B\\H*(w{E^0))) + (w{E^0)) + l) [ {Bu(t),F(u(t))) x dt. (18) 

Jo 

Proof. To prove this inequality, we use the isometric representation of B in the space L 2 (fl,p). 
Denoting / = U~ 1 u, using that U~ 1 BUf = bf , p(f) = U~ 1 F(U /), we have, for instance, 
||S 1 / 2 it||^- = (/, 6 /)l 2 (n,/^)j an d hence it suffices to prove that 


w(E*(0)) / ( bf 2 + bp(f) 2 )dpdt 


< 


T [ bdpH*(w(E^0))) + (w{E^(0)) + l) [ [ bfp(f)dpdt. (19) 

J n Jo Jn 


Indeed, this implies (18) (note that bdp ^ ||S|| by the spectral theorem). 

Let us prove (19). First of all, for every t G [0,T] we set U* = {x G fl | \f(t,x)\ ^ £o}- If 
b = 0 on fl* then the forthcoming integrals (see in particular the left-hand side of ( 20 ) are zero 
and there is nothing to prove; hence, without loss of generality we assume that b is nontrivial on 
Q*. Using ( 6 ), we choose eo > 0 small enough such that -pip{f ) 2 ^ Sq almost everywhere in S2‘, 
and therefore we have 


Sn\ bd P 




—p(f) 2 bdp G [0,sg 


Using the Jensen inequality with the measure bdp , and using the fact that H(x) = \fxg{\fx), we 
get 


H 


In* bd V 


I 4 p(f) 2bd p] < j \ , / —\p(f)\g(—\p(f)\]bdp. 

Jan <4 / l 0 , bdp. J n t c -2 \c 2 J 


Using ( 6 ), we have |p(/)(x)| ^ C 2 g 1 (|/(a;)|) for almost every x G f2‘, and since g is increasing, 
we get that g I d-\p(f)\ ) < |/| almost everywhere in f}*. Since fp(f) ^ 0 by (5), we get 


H 


\T~rr [ ~aP(f) 2bd A ^ T~ 

\ in* b d d Jn\ c 2 J Jo* 


1 1 


b dp C2 

“ -“i 

Since FJ is increasing, it follows that 

j bp(f) 2 dp ^ c 2 [ bdp H^ 1 
J n; J n* 


( bifMDidp^T^-- f 

n* b dp a J a 


b fp(f) d P- 


_J_I f 

ju\bdp c 2 J n 


b fp{f) d P , 
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and therefore 
r T 


/o 


w{E+(0)) [ bp(f) 2 dp dt ^ [ wiE'fri 0))c 2 [ bdp H 1 ( , \ — f bf p(f) dp] dt. 

Jn\ Jo Jn\ yjfjtbdpua Jn J 


Thanks to the Young inequality AB ^ H(A) + H*(B) (where H* is the convex conjugate), we 
infer that 


7 (^( 0 )) f bp{f) 2 dp dt ^ f c 2 [ bfp(f) dp + f cl f bdp H*(w(E^(Q))) dt 
J Q\ J 0 J Q, Jo Jn\ 


( 20 ) 


< C 2 / / bfp(f) dpdt + c 2 T / bdp H*{w(E cj> { 0))). 

Jo Jn Jn 


Besides, in 12 \ 12^, using ( 6 ) we have \p{f)\ < |/|. Using (5), it follows that 

\2 


f w{E <j> { 0)) [ bp(f) 2 dpdt< f w(E <j) ( 0)) [ b\f\\p(f)\dpdt 
0 Jn\n\ Jo Jn\n\ 

< J w{E (j> { 0 )) [ bf p(f) dp dt. 


( 21 ) 


in 


From (20) and (21), we infer that 


[ w(E^0)) f bp(f) 2 dpdt<{w(E4 0)) + l) f [ bfp(f) dpdt+ T [ bdp H*{w(E^( 0))). 

Jo Jn Jo Jn Jn 

Let us now proceed in a similar way in order to estimate the term w(E c / > ( 0)) J^bf 2 dpdt. 

We set rf = iL _1 and £1 = min(so, g(r\)) < 1. For every t 6 [0,T], we define 121 = 

{x 6 12 | |/(2, x)| ^ £ 1 }. As before, without loss of generality we assume that b is nontrivial on fl^. 
From (6), we have cig(\f\) ^ \p(f)\ in 12|- By construction, we have 

f ! , [ f 2 bdp£[0,sl}. 

Jn* 6< Wn* 

Using the Jensen inequality as previously, and using (6) and (5), we infer that 

1 r . o, , \ . 1 


H 1 -i—et / ^ < 717 / l/llff(/)I^M 

Jq* bdp Jnt 2 J JntbdpJn* 




1 


cif^ bdp J n * 


&l/l|p(/)|dju < 


1 


ci Ini bd pJn 


bfp(f)dp 


Since H is increasing, and integrating in time, we get 
r T p pT 


Jo 


w(E c j > { 0)) [ f 2 bdpdt ^ [ w(E<j,{ 0)) [ bdp H 

J Jo J cii 


Cl S a t bdp Jn 


bfp(f) dp dt 


It then follows from the Young inequality that 

r T 


/ 0 


w(E ( j,{ 0)) f f 2 bdpdt f bdpH*{w{E${ 0))) H- [ [ bfp(f) dpdt 

Jnt. Jn ci Jo Jn 


The estimate in 12 \ 12^ is obtained in a similar way. 
The lemma is proved. 


□ 
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Fourth step. End of the proof. 

Lemma 4. We have 

E U (T) < E u ( 0) ^1 - p T L~ 
for some positive sufficiently small constant px ■ 


-i (E u { 0) 


( 22 ) 


Proof Using successively the observability inequality (9), the estimate (16) of Lemma 2 and the 
estimate (14) of Lemma 1, we first get that 

20^(0) < [ \\B 1,2 (j>{t)\\ 2 xdt^k T f \\B^ 2 z(t)\\ 2 x dt 

Jo Jo 

< 2/ct £ (||B 1/ Vt)|H- + II B^Ffumix) dt. 

Multiplying this inequality by the constant w(E t f,( 0)), it follows from the estimate (18) of Lemma 
3 that 

C T w(E^0))E < t > (0)<kTT\\B\\H*(w{E <l> (0))) + kT{w(E <l> (0)) + l) [ {Bu(t),F(u(t))) x dt, (23) 

Jo 

From (7) and (17), we have L(w(s)) = -| = H for every s £ [0,/3sq), and hence (3H*(w(s)) = 

sw(s). We choose ft large enough such that E^,( 0) < fts ^, and thus in particular we get 


H*(w(E^ 0))) = 


w{Et( 0 ))E 40 ) 

13 

E <t>( o) 


(24) 


Besides, we also choose ft large enough such that /3 ^ l(H'(s 2 )) > an< ^ s i nce L '■ [0, +oo) —> [0, Sg) is 
continuous and increasing, it follows that 


w(E^0)) < H'(sl) < 1. 
Finally, we get from (23), (24) and (25) that 
k T T\\B\\\ 


(25) 


C T - 


13 J 


w(E^,(0))E^(0) < k T / (Bu(t),F(u(t))) x dt, 


We choose ft large enough such that Ct — kTT ^ B ^ > It follows that 


Ptw(E^(0))E^(0) < / ( Bu(t ), F(u{t))) x dt = £ u (0) - E U (T), 


(26) 


for some positive sufficiently small constant px, and since £^(0) = E u ( 0), using (17), we have 
obtained that 

E U (T) < E u { 0) (l - pxL- 1 

□ 


as expected. 

By translation invariance, we get from (22) that 

E u {{k + 1)T) < E u {kT) f 1 - Pt L 


E, 




13 JJ ’ 
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for every integer k. Setting Ek = and M(x) = xL~ x {x) for every x £ [0, s§], we have 

Ek+i — Ek + pxM^Ek) ^ 0, for every integer k. From these inequalities, we deduce by routine 
arguments (similar to those in [9], and thus not reproduced here) that, setting K r (r) = M( y ) > 
we have 


M(E p 


^ — min 
Pt ^e{o,...,p} 


K r 1 (p T (p-i)) 
£ +1 


from which it follows that 


E -<‘>« e ™- 1 („<&., G*" 


<0 TL 



for sufficiently large t, with ijj defined by (8). The precise constants in the estimate (10) follow 
from the choice of /? (large enough) above. Theorem 1 is proved. 


Proof of Corollary 1. Since </(0) ^ 0 and g is strictly increasing, there exist ci > 0 and C 2 > 0 
such that 

ci|s| ^ \g(s)\ < c 2 |s| V |s| < 1. (27) 

To overcome the difficulty of dealing with non strictly convex functions, let us consider for e > 0 
small enough, the function g e (s) = s 1+£ for |s| ^ 1. According to (27), we have 

ci 9e(\fi x )\) < \p(f)i x )\ < c 2 g~ 1 {\f(x)\) for almost every x eCl such that |/(x)| < 1, 


Denote respectively by L e and w e the functions defined from g = g e respectively by (7) and (17). 
Straightforward computations lead to 


L e{r) 


£ 

2 + £ 



L-\r) 


2 + e ^(2 + e)r^ 2 


for |r| ^ 1 and it follows that 


w E {s) = L e 1 



2 + e /(2 + e ) S y /2 

2 V Pe ) 


for s small enough. Hence, the family (w e ) e> o converges pointwisely on (0,/3sg) to 1 as e tends to 

0 . 

Following the proof of Theorem 1, one shows that the estimate (26) is verified with w E as weight 
function, in other words that 


p T w e (E <t ,(0))E < i > (0) ^ E u ( 0) - E U (T), 

for the same positive constant pt £ (0,1) as the one introduced in (26). Letting e go to 0 and 
using that E^(Q) = E u { 0) leads to 

0 <E U (T) < {1- Pt )E u (0). 

It is standard to derive from such an estimate the exponential decay of the energy E u (see e.g. 
[21]). For that purpose, let us reproduce the previous reasoning on each interval [jT. {j + 1)T) 
with j Js 1. Using an induction argument, it follows that 

0 < E u (jT) < (1 - Pt ) j E u ( 0) = e~ uj E u ( 0), 

where ui = — ln((l — pr)) > 0. Then, for every t > T, there exists a unique j £ N such that 
t £ [jT, (j + 1)T). Since j > T — 1, and E is nonincreasing, we obtain 

0 < E u (t) < E u {jT) < e~ uj E u { 0) = -^— e -^ T E u { 0), 

1 ~ Pt 

whence the expected result. 
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3 Discretization issues: uniform decay results 

3.1 Semi-discretization in space 

In this section, we introduce and analyze a general space semi-discrete version of (1). Our main 
objective is to prove a theorem similar to Theorem 1, but in this semi-discrete setting, with 
estimates that are uniform with respect to the mesh parameter. In order to ensure uniformity, we 
add an extra numerical viscosity term in our discretization scheme (as in [40]), and we establish 
uniform decay rates estimates as those proved in Theorem 1. 

Let Ax > 0 be a space discretization parameter, standing for the size of the mesh. We assume 
that 0 < Ax < Ax 0 , for some fixed Ax 0 > 0. We now introduce an appropriate setting for 
semi-discretizations in space, following [25, 27] (see also [14]). Let (Xai)o<Ai<Ai 0 be a family of 
finite-dimensional vector spaces. Without loss of generality, we identify X Ax with ]FI Ar ( Ax ), where 
N{Ax) € IN. More precisely, if {xi} i^i^n(Ax) is a family of mesh nodes and a is a mapping 
from X to 1R, then the vector u A x £ Xai is expected to be an approximation of the vector 
('u(xd)l^i^iV(Ax) • 

Hereafter, the notations < and ~, already used in the continuous setting, keep the same meaning 
as before, with the additional requirement that they also mean that the involved constants are 
uniform as well with respect to Ax. 

Let us recall what is the Hilbert space usually denoted by X 1 / 2 (see, e.g., [17]). Let /3 be 
a real number belonging to the resolvent of A. Then the space Xj / 2 is the image of X under 
(/3idjf — A) -1 / 2 , and it is endowed with the norm ||u||x 1/2 = ||— AY^uWx- For instance, 
when A 1 ' 2 is well defined, we have X\/ 2 = D^A 1 ' 2 ). We set X_i/ 2 = X[/ 2 , where the dual is 
taken with respect to the pivot space X. 

We assume that, for every Ax £ (0, Ax 0 ), there exist linear mappings Pax '■ X_i/ 2 —•► X Ax 
and P Ax : X Ax —^ X \/2 such that PaxPax = idx Ax ■ We assume that the scheme is convergent, 
that is, ||(I — P A xPax)u\\x —> 0 as Ax —>• 0, for every u £ X. Here, we have implicitly used 
the canonical injections D(A) °->- X\/ 2 X X_i/ 2 (see [17]). Additionally, we assume that 
Pax = P A x- Note that, at this step, these assumptions are very general and hold for most of 
numerical schemes. 

For every Ax £ (0, Ax 0 ), the vector space X Ax is endowed with the Euclidean norm || ||ax 
defined by ||max||ax = ||-PaxWAx||.y, for u A x £ X A x- The corresponding scalar product is denoted 
by (*, *)ax- Note that [I-PaxIIhWa*^) = 1 and that, by the Uniform Boundedness Principle, 
\\Pax\\l(x,x^) < L 

For every Ax £ (0,Ax o ), we define the approximation operators Aax '■ Xax —> Xax of 
A, and BAx '■ X Ax X A x, by A A x = PaxAPax (where, in this formula, we have implicitly 
used the canonical extension of the operator A : X\/ 2 —> X_i/ 2 ) and Bax = PaxBPax- Since 
Pax = FA, note that A Ax is identified with a skew-symmetric matrix, and B Ax is identified with 
a symmetric nonnegative matrix. 2 Finally, we define the (nonlinear) mapping Fax '■ X A x —> X A x 
by Fa x (uax) = PaxF(PaxUax), for every u A x £ X Ax - Note that, by construction, B Ax is 
uniformly bounded with respect to Ax, and Fax is Lipschitz continuous on bounded subsets of 
Xax, uniformly with respect to Ax. 

We consider the space semi-discrete approximation of (1) given by 

u'ax^) + A AxUAx{t) + B Ax FAx(uAx(,t)) + (Ax) CT Vax^Ax W = 0. (28) 

The additional term (Ax) a VAxUAx{t), with a > 0, is a numerical viscosity term whose role is 
crucial in order to establish decay estimates that are uniform with respect to Ax. The role and the 

We have (BaxUAx ? Ax — — {BPaxUAx.^ Pax.u aA x x- 0. 
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design of such a term will be discussed further. We only assume, throughout, that Vax ■ —> 

Xai is a positive selfadjoint operator. 

Defining the energy of a solution u A x of (28) by 

,{P) = ~ll u Ax(OII Ax; (29) 

we have, as long as the solution is well defined, 

E'u Alc (t) = ~{uAx(t) : B Ax F Ax (uAx{t)))Ax - (&xY\\(VAx) 1/2 UAx{t)\\ 2 Ax . (30) 

We are going to perform an analysis similar to the one done in Section 2, but in the space semi- 
discrete setting, with the objective of deriving sharp decay estimates for (28), which are uniform 
with respect to Ax. 

3.1.1 Main result 

Assumptions and notations. First of all, we assume that 

{uAx,BaxFax(uAx))ax + {Ax) a \\(VAx) 1/2 UAx\\Ax > 0, (31) 

for every u A x £ X A x- Using (30), this assumption ensures that the energy E uz ^ x (t) defined by 
(29) is nonincreasing. 

Moreover, we assume that there exists so > 0 such that 

sup {AxYWV^W^x^ < +oo. (32) 

Ax6(0,so] 

This assumption will be commented in Remark 3. 

We keep all notations and assumptions done in Section 2.1. In order to discretize the mapping 
p : L 2 (D, fj,) —> L 2 (Q,p) defined by p(f) = U -1 F(U /), we use as well the approximation spaces 
Xax, as follows. We first map the functional setting of X to L 2 (fl,p) by using the isometry 
[/ -1 : X —> L 2 (tt,p). In Section 2.1, we have defined the operator A = U~ 1 AU on H = L 2 (n,p), 
of domain Hi = D{A) = U~ 1 D(A). Accordingly, we set Hi /2 = U~ 1 X 1 / 2 , and we define %_ i / 2 as 
the dual of H \/2 with respect to the pivot space H. We have H_i/ 2 = U~ 1 X_ij 2 , where we keep 
the same notation U to designate the canonical (isometric) extension U : H_\/ 2 X_i/ 2 . Now, 

for every Ax £ (0,Axq), the linear mappings P Ax U : H_ i / 2 —> X Ax and U~ 1 Pa x ■ X Ax —> H\/ 2 
give a space discretization of the Hilbert space H = L 2 (fl, p) on the finite-dimensional spaces X A x- 
For every u A x £ Xa x , we set 

PAxi'U'Ax') = PAxU p{U PAx^Ax) = PaxF(Pax'UAx) — Fax Yax ): 

so that, finally, we have pAx = Fax- 

Besides, for every / G L 2 (D,/z), we set 

PAx(f) = U^PaxPaxUpU) = U^pAxPAxFiUf). (33) 

By definition, we have Pax(U~ 1 PaxUax) = U~ l P A x PAx{uax), for every uax £ X Ax - The map¬ 
ping pAx is the mapping p filtered by the “sampling operator" U^ 1 PaxPaxU = (PaxU)*PaxU. 
By assumption, the latter operator converges pointwise to the identity as Ax —> 0, and in many 
numerical schemes it corresponds to take sampled values of a given function /. 

We have pAa:(0) = 0, but (5) and ( 6 ) are not necessarily satisfied, with p replaced with pAx- 
In the sequel, setting fAx = U~ 1 PaxUax, we assume that 

fAx P Ax (fAx) > 0, (34) 
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and that 


ci g(\fAx{x)\) < \pAx(fAx)(x)\ s? c 2 g '(|/a*(®)|) for a.e. x G fl such that \fAx(x)\ ^ 1, 
ci |/a* 0)| < |pA a; (/Ax)(a;)| < c 2 |/axO)I for a - e - ® G such that |/axO)I ^ 1) 

for every uax £ Xah for every Ax € ( 0 , Axo). 

Note that the additional assumptions (34) and (35) are valid for many classical numerical 
schemes, such as finite differences, finite elements, and in more general, for any method based on 
Lagrange interpolation, in which inequalities or sign conditions are preserved under sampling. But 
for instance this assumption may fail for spectral methods (global polynomial approximation) in 
which sign conditions may not be preserved at the nodes of the scheme. The same remark applies 
to (31). 

Note also that assuming (34) and (35) is weaker than assuming (5) and ( 6 ) with p replaced 
with pax, for every f £ L 2 (D,p). Indeed, the inequalities (34) and (35) are required to hold only 
at the nodes of the numerical scheme. 

The main result of this section is the following. 

Theorem 2. In addition to the above assumptions, we assume that there exist T > 0, a > 0 and 
Ct > 0 such that 

CtE^M < £ (lIB^AxWIlL + (Ax) CT ||Vi / ^AxWllL) dt, ( 36 ) 

for every solution of (f>'^ x {t) + AAxf>Ax(t) = 0 (uniform observability inequality with viscosity for 
the space semi-discretized linear conservative equation). 

Then, the solutions of (28), with values in Xax, are well defined on [0,+oo), and the energy 
of any solution satisfies 

Eu^it) ^ Tmax (.hi E u ^ x (0)) L 

for every t > 0 and Ax 6 (0, so], with 71 ~ || 5||/72 and 72 — Ct/T{T 2 \\B 1 / 2 \\ 4l + 1). Moreover, 
under ( 11 ), we have the simplified decay rate 

Eu±.(t) < T max( 7 i, E u&x (0)) {H'y 1 (j) , 

for every t > 0 and Ax 6 (0,so], for some positive constant 73 ~ 1 . 

Remark 2 (Comments on the uniform space senri-discrete observability inequality (36)). The 
main assumption above is the uniform observability inequality (36), which is not easy to obtain in 
general. There are not so many general results in the existing literature, providing such uniform 
estimates. 

First of all, in the absence of a viscosity term, the observability inequality (36) fails to be uniform 
in general (see [)5] and references therein), in the sense that the largest constant appearing at the 
left-hand side of (36), depending on Ax in general, tends to 0 as Ax —> 0. This phenomenon, 
which is by now well known, is due to the fact that the discretization creates spurious highfrequency 
oscillations that cause a vanishing speed of highfrequency wave packets. We refer to [)5] for the 
detailed description of this lack of uniformity, in particular for ID wave equations with boundary 
observation (for which a simple computation yields non-uniformity). 

A counterexample to uniformity is provided in [31] for the ID Dirichlet wave equation with in¬ 
ternal observation (over a subset ui), semi-discretized in space by finite differences: it is proved that, 
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for every solution of the conservative system 4>'& x {t) — Aai^a xif) = 0 , there holds Cai-E^, ( 0 ) ^ 

Jo XuM'AxitfWAx where the largest positive constant Cax for which this inequality is valid sat¬ 
isfies C Ax —» 0 as Ax —> 0. The proof of this lack of uniformity combines the following facts: 
using gaussian beams, it is shown that, along every bicharacteristic ray, there exists a solution 
of the wave equation whose energy is localized along this ray; the velocity of highfrequency wave 
packets for the discrete model tends to 0 as Ax tends to 0. Then, for every T > 0, for Ax > 0 
small enough, there exist initial data whose corresponding solution is concentrated along a ray that 
does not reach the observed region oj within time T. 

Note that the uniform observability inequality (36) holds for all solutions of the linear conser¬ 
vative equation f Al (t) + Aax4 > Ax(() = 0, if and only if the solutions of z' Ax (t) + AAxZAx(t) + 
BAxZAxit ) = 0 are exponentially decaying, with a uniform exponential rate (see [19], see also 
Lemma 6 and the end of Section 3.1.2, from which this claim follows by an easy adaptation). 

Many possible remedies to the lack of uniformity have been proposed in [45], among which the 
filtering of highfrequencies, the use of multigrids, or the use of appropriate viscosity terms. Here, 
since we are in a nonlinear context, we focus on the use of viscosity, that we find more appropriate. 3 
The role of the numerical viscosity term (which vanishes as the mesh size tends to zero) is to damp 
out the highfrequency numerical spurious oscillations that appear in the semi-discrete setting. 

A typical choice of the viscosity operator, proposed in [19], is 

Vax = A Ax Aaxi 

(symmetric positive definite square root of AaxI if it is not positive definite, add some eid.Y Aa J, with 
a chosen such that (31) is satisfied. Some variants of the viscosity term are possible. Unfortunately, 
it is not proved in the existing literature that such a general viscosity term is systematically sufficient 
in order to recover the desired uniform properties (in contrast to time discretizations, see further). 
As discussed in [19], the main difficulty consists in establishing the uniform observability inequality 
(36) with viscosity (see [45] for a thorough discussion on this issue). 

There are quite few results in the literature where one can find such uniform stability results, 
for some particular classes of equations. 

One of them concerns the wave equation, studied in ID (and in a 2D square) in [40], with an 
internal damping (see also the generalization to any regular 2D domain in [35]), semi-discretized 
in space by finite differences. Using discrete multipliers, the authors prove that the solutions of the 
semi-discretized locally damped wave equation with viscosity 

y'kxtt) - Aa*2/a x {t) + a A xVAxty) ~ (Ax) 2 A a xV'axW = 0, 

decay exponentially to 0 ; uniformly with respect to Ax. Here, the viscosity term is — (Ax) 2 AaxVax’ 
and the term oai stands for the discretization of a localized damping. With our notations, we have 


UAx{t) - few) ’ 


( 0 -I N \ 

\-Aax 0 ) ’ 




(° «!)• v -=(o -D 


with cr = 2, and with the usual finite-difference discretization of the Laplacian, given for 

3 Note that, in the linear context, filtering and adding a viscosity term are equivalent, as it follows from [19, 
Corollary 3.8 and Remark 3.9]. 
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instance in ID by 


Aax — 


{Ax) 2 


(-2 1 
1 
o 

V o ... 


°\ 


0 1 - 2 / 


In that case, the viscosity operator is not positive definite, however the uniform observability in¬ 
equality (36) follows from the proof of [19, Theorem 1.1]. 

In [18, Theorem 7.1], the author studies a class of second-order equations, with a self-adjoint 
positive operator, semi-discretized in space by means of finite elements on general meshes, not 
necessarily regular. This result is generalized in [34] with a weaker numerical viscosity term. We 
refer also to [1] for a similar study under appropriate spectral gap conditions, noting that, in that 
paper, results are also provided on uniform polynomial stability. 

Finally, in [38] and [39], the authors use a similar viscosity operator for the plate equation 
and even for a more general class of second-order evolution equations (under appropriate spectral 
gap assumptions, and thus essentially in ID). They do not explicitly prove a uniform observability 
inequality but they establish a close result based on a frequential characterization of the dissipation 
introduced in [29]. 


Remark 3 (Comments on the assumption (32) on the numerical viscosity.). The assumption (32) 
is satisfied in all cases mentioned in Remark 2 (see the corresponding references): in [40], some 
explanations about the choice of the viscosity operator and its properties (including (32) ) are given 
after Theorem 1.1; in [19], (32) corresponds to the third assumption on the viscosity operator in 
Theorem 3.7; in [35], a finite difference approximation scheme is used, with a = 2 in the viscosity 
(it is easy to see that (32) is indeed satisfied); finally, in [38, 39], this assumption follows from 
their particular choice of the numerical viscosity operator. 


3.1.2 Proof of Theorem 2 

We follow the lines of the proof of Theorem 1. The only difference is in handling the viscosity term. 
Hereafter, we provide the main steps and we only give details when there are some differences with 
respect to the continuous case. 

Note that the global well-posedness of the solutions follows from usual a priori arguments (using 
energy estimates), as in the continuous setting. Uniqueness follows as well from the assumption 
that F is locally Lipschitz on bounded sets. 

First step. Comparison of the nonlinear equation (28) with the linear damped model. 

We first compare the nonlinear equation (28) with the linear damped system 

z Ax(t) + ^Ai*Ai(f) + BAxZAx{t) + {Ax)° VAx z Ax{t) = 0. (37) 

Lemma 5. For every solution uax{') of (28), the solution zA x(') of (37) such that zax( 0) = 
UAi(O) satisfies 

£ (K /2 w*)iiL + (Axnivi 7 x 2 wi)iiL) dt 

< 2 £ (\\BTxUAx{t)\\l x + WB^FA^UAximlx + 2(Axn|V^Vx(t)||L) dt. 
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Proof. Setting 0 a x (t) = u& x (t) - z& x (t), we have 


(^Ax(^) "h -^-Ax^Axif) H” ^Ax^Ax(^Ax(^)) ^Ax^Ax(^) H” (A#) ^Ax^AxijP) ^Axifpi) Ax — 0- 

Denoting £^ Ax (t) = |||V>Ax(t)|| Ax , it follows that 

^| A ,W + ll-®Ax' ZAa: (^)ll Ax + (A*) II^Ax^AxIIa® = ~( u Ax{t), B/\ x F^ x (uAx{t))) Ax 

+ (-^Ax-^Ax^AxW). S&j^AxWJa! + (B^UAxit), B^ Z A x(t)) Ax ■ 

Using (31), we have (u Ax (f), B Ax .F Ax (u Ax (f))} Ax > —(Aa;) 0- 1| V^ 2 M Ax (f)||^ x , and hence 

^0 A .W + II^AxWIIL + (AxJHIV^VaxIIL < l|5i / ^Ax(« A xW)|| A x|| J Bi / A A x(i)l|Ax 

+ \\BiluAx{t)\\Ax\\Bi 2 x ZAx{t)\W + (A^HIViCVxWIlL- 

2 >2 -i 

Thanks to the Young inequality a& ^ §e + with 9 = we get 

^ A .(t) + II^aaxWIIL + (A^HIvi^AxIlL 

< \\\ B TxZAx(t)W\x + \\BT x FAx{UAx(t))\\\x + II^VxWIlL + (A*)" || V^Vx («) ||L. 
and thus, 

£* A .W + \\\BTx^x{t)Wlx + (AxJHIVi^AxIlL 

< WB^FAxMmix + ws^uAxmix + (A^-nv^VxWiiL- 

Integrating in time, noting that E^ &x (0) = 0, we infer that 

E^(T) + \ £ (WB^ZAxmlx + (A*) ff ||V^ 2 z Al ||L) dt 

< £ (WB^FAxiuAxmlx + II^xVxWIIL + 2(Ax)H|vi / x \ A x(t)||L) dt. 

Since E^^ x (T) ^ 0, the conclusion follows. □ 

Second step. Comparison of the linear damped equation with the conservative linear equation. 
We consider the space semi-discretized conservative linear system 

0 Ax (t) + Aax4>A x{t) = 0. (38) 

Lemma 6. Assume that (32) holds true. Then, for every solution zax(‘) of (37), the solution 
4 >Ax °f (38) such that <(> Ax (0) = -z Ax (0) satisfies 

[ (ii^x^axWiil + (Az)ivi^ Ax (t)iiL) dt 

< hr £ (||<AAxWIIL + (AxYWV^ZAxmlx) dt, 
with k T = 1+T 2 +T 2 \\B 1 / 2 \\ 2 +T 2 \\B 1 / 2 \\ a . 
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Proof. Setting 9 Ax (t) = <t> Ax {t) - z Ax (t) and E g ^ x (t) = ^\\9 Ax (t)\\ 2 Ax , we have 

(0 Ax (t) + A Ax 9 Ax - B Ax z Ax (t) - ( Ax) a V Ax z Ax (t ), 0 Ax [t)) Ax = 0, 

and therefore 

E'e^ x (t) = ( B Ax z Ax (t),9 Ax (t)) Ax + {Ax) a (V Ax z Ax {t),9 Ax (t)) Ax . 

Integrating a first time over [0,i] and a second time over [0,T], and noting that Eg^ x ( 0) = 0, we 
get 

\\0 Ax {t)\\ Ax dt = J (T-t)({B Ax z Ax {t),9 Ax (t)) Ax +(Ax) a (V Ax z Ax {t),9 Ax (t)) Ax ) dt. 

Applying the Young inequality ab ^ + A for some v > 0 to both terms at the right-hand side 

of the above equality, we get 

pT r T 


'0 


pAx(t)\\ Ax dt ^ T 2 v [ \\B Ax z Ax (t)W 2 Ax dt+ - [ \\9 Ax {t)W 2 Ax dt 
Jo v Jo 

+ T 2 HAx)° [ T \\V Ax z Ax (t)\\ 2 Ax dt + r^ x 9 Axm 2 Ax d, (39) 

Jo v Jo 


Using moreover that 

r T 


(nwt)iiL + (^r\\^ x e Ax mA X ) *< m Ax J^ u«w*)iil^ 

with M Ax = 1 + (AxYWV^lW^x^y one gets for all v > M Ax 

MaA / ' T (ii0A,wiiL + (A^nivi / >A,wiiL) dt 

^T 2 uM Ax [ \\B Ax z Ax (t)\\ 2 Ax dt + T 2 v(Ax) a M Ax f \\V^z Ax (t)\\ 2 A x dt. 
Jo Jo 


1 - 


/ 0 

Let us choose v = 2Ma x , the last inequality becomes 

1 rT 


r 0 


(llflAx(t)HL + (Axr\\V^9 Ax ml x ) dt 

< 2 T 2 M 2 Ax f f \\B Ax z Ax (t)\\ 2 Ax dt + (A x y £ IIV^zaxWIIL dt ) ■ 


As a result, there holds 

u T 

2 jo 


(llflAx(t)HL + (A*) ff ||V^0Ax(t)||L) dt 

< r 2 ||S 1/2 || 2 f T \\B^z Ax ml x dt + T 2 (AxY [ T ||V^Ax(t)||L dt, 

Jo Jo 
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where, to obtain the latter inequality, we have used that ||i?Ax|| ||-B||. Writing <f>Ax = ^ai + ^Ah 

we have 



{\\B^<j)Ax(t)\\A x dt + (A^^IIVa^AxWIIL) dt 

(\\B^0A X mi x dt+{Axnv^e^mi x ) 


^ 2 


dt 


+ 2 



[\\B)llz Ax {t)\\\ x dt + (AaOHIVi^AxWIlL;) dt, 


and the lemma follows. 


□ 


Third step. Nonlinear energy estimate. 

We define w by (17), as before, with /3 > 0 to be chosen later. 

Lemma 7. For every solution u A x of (37), we have 


[ (K / *V*(i)||L + \\BT x F Ax {u Ax m\ x + (Aa:) cr ||V^UAx(i)|| Ax) dt 

<T\\B\\H*(w(E^m 

+ (wiE^iO)) + 1) J ^{B A xUAx{t),F Ax (uAx(t))}Ax + {^x) a \\V^UAx{t)\\ 2 Ax ^ dt. 

Proof. We set f A x(t) = U~ 1 P A xU A x(t). Using the definitions of the discretized operators and the 
isometry U, we have 


|| ^ Ax 'U'Ax [| Ax — {Bax^Ax, aAx) Ax — PAx^Ax•> PAx^Ax)X 

~ {U BUU Pax'OAxtU Pax^Ax) L 2 (Q,/i) = {bf A x, /*Ax)l 2 (Q,/a) = / ^/aj; 

Jn 

and, using (33), 

|| B Ax ^Ax (^Ax) || Ax = Ax F Ax [a Ax) > PAx (^ Ax)) Ax = (.-BPaxEA x (^Ax)j PAx EAx (^Ai)) X 

= {U^ 1 BUU~ 1 PaxPAxF(Pa x UAx), U” 1 PaxPAxF(PaxUAx)) L 2 (fi,/i) 

= {bU -1 PaxPaxF(UU~ 1 PaxUAx), U~ 1 pAxPAxF(UU~ 1 pAxUAx))L 2 (n,iJ.) 
= {bpAx(fAx), PAx(/Ax))l/ 2 (Q,/A) 

= [ b(pAx{fAx)) 2 d/I, 

Jn 

and, with a similar computation, 

{Bax'O'Ax) EAx {f Ax)) Ax — j bfAxPAx^fAx^dp. 

Jn 

Then, to prove the lemma, it suffices to prove that 

J w(-M°))^ (bfAx + b(pAx{fAx)) 2 ^j dpdt 

<T f bdp FJ* {^(E^if)))) + (w(i^(0)) + 1) [ [ 6/axPAx(/ax) dp dt. 

J Q JO J Q 
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Indeed, the viscosity terms at the left-hand and right-hand sides of the desired inequality play 
no role in that lemma. Then, the proof is completely similar to the one of Lemma 3, using in 
particular the assumptions (34) and (35). We skip it. □ 


Fourth step. End of the proof. 
Lemma 8. We have 


E U&X (T) < E u&x (0 ) (l -p T L- 1 


for some positive constant pt < 1. 

Proof. Recall that £^ Aa , (0) = E u&x (0) and that 

^( 0 )-^.(r)= f ((B Ax u Ax (t),F Ax (u Ax (t))) Ax + (Axr\\V 1 ^u Ax (t)\\l x )dt. ( 40 ) 

Jo 

Using successively the observability inequality (36) and the estimates of Lemma 6, of Lemma 5 
and of Lemma 7, we get 


2C T w{E^M)E^M 

(lIs^AxWllL + CAxrilvi^AxWIlL) dt 
<k T w(E^ x ( 0))^ T (||B^x(t)||L + (A*ri|V^AxW||L) dt 

< 2fc T w(F 0Ax (O)) J {\\B^u Ax (t)\\ 2 Ax + \\B][lF Ax {u Ax {t))\\ 2 Ax + (A:r) CT || V^ x u Ax (t)\\ L,) dt 
<k T T\\B\\H*(w(E^ x m 

+ {w(E < / )&x ( 0)) + l) J [(F Ax (u Ax (t)),B Ax u Ax (t)) Ax + (Ax)' T ||V^ c 2 MA a: (i)||Ax) dt. 


Using (40), we infer that 


2C t w{E^M)E^M < k T T\\B\\H*(w(E ( f >&x (0))) + («,(£* A .(0)) + l)(F„ Ax (0) - E u&x (T)). 
We conclude similarly as in the proof of Lemma 4. □ 


The end of the proof of Theorem 2 follows the same lines as in the continuous case. 


3.2 Semi-discretization in time 


In this section, we analyze a time semi-discrete version of (1), with the objective of establishing 
uniform decay estimates. To this aim, following [19], we add a suitable viscosity term in an implicit 
midpoint numerical scheme. 

Given a solution u of (1), for any At > 0, we denote by u k the approximation of u at time 
tk = kAt with k G IN. We consider the implicit midpoint time discretization of (1) given by 


u k +i _ u k 

At + 

~k+l _ u k+1 

At 

K u°=u( 0). 


u k + u k+1 


+ BF 


u k + u k+1 


= 0 , 


= V At u k+ \ 


(41) 
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The second equation in (41) is a viscosity term (see Remark 4 further for comments on this term 
and on the choice of the midpoint rule). We only assume, throughout, that Va t ■ X —> X is a 
positive selfadjoint operator. 

Written in an expansive way, (41) gives 


u k+1 - u k 
A t 


+ A 


./c +1 


BF 


+ (id A - + A tV At )u k+1 


+ V At u k+1 + ^AV At u k+1 = 0. 


We define the energy of a solution (u k )k^TN of (41) as the sequence (E u k)ke^s given by 

1 , 


E„h = - m 


k || 2 


IX- 


As long as the solution is well defined, using the first equation in (41), we have 

,F 


Efik+i — E u k = —At l\B 


u k + u k+1 


u k + u k+1 


x 


and using the second equation in (41), we get that 

E&+* =E u k +1 +At\\(V At ) 1 '\ k+1 \\% + ^\\V At u k+l fx- 

We infer from these two relations that 


(42) 


- /\1~' 




u k + u k+1 \ \ 

2 / lx 

- At II (VaO 17 V+ 1 III - ^ II V At « fe+1 III, 


(43) 


for every integer k. Note that, thanks to (4), we have E u k+i — E u k ^ 0, and therefore the energy 
defined by (42) decays. We next perform an analysis similar to the one done in Section 2, but in 
the time semi-discrete setting, with the objective of deriving sharp decay estimates for (28), which 
are uniform with respect to At and Ax. 


3.2.1 Main result 

We keep all notations and assumptions done in Section 2.1. 

Hereafter, the notations < and already used in the previous sections, keep the same meaning 
as before, with the additional requirement that they also mean that the involved constants are 
uniform as well with respect to At. 

The main result of this section is the following. 

Theorem 3. In addition to the above assumptions, we assume that there exist T > 0 and Ct > 0 
such that, setting N = [T/At] (integer part), we have 


N -1 


CtUYx < a tJ2 


k=0 


B x ' 2 


A*:+l 


X 

N -1 


N -1 


+ Ai £ ||(V At ) 1/2 ^ +1 Hx + (At) 2 £ ||V A ^ fc+1 ||!, (44) 


k—0 


k—0 
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for every solution of 


> k 




U k + 1 - /A fc ^ ^ \ ^ 


At 

A^+l _ ✓/»&+! 


( 45 ) 


Ai 


= V A t# 


k+1 


(uniform observability inequality with viscosity for the time semi-discretized linear conservative 
equation with viscosity). 

Then, the solutions of (41) are well defined on [0, +oo) and, the energy of any solution satisfies 


for every integer k, with 72 ~ Ct/T( l + e 2T ^ B ^ max(l, T||B||)) and 71 ~ ||U||/72. Moreover, under 
(11), we have the simplified decay rate 


E uk < T max( 7l , E u o ) ( H 0" 1 , 

for every integer k, for some positive constant 73 ~ 1. 

Remark 4. As in Remark 2, we insist on the crucial role of the viscosity. 

In the absence of a viscosity term, the decay is not uniform in general: see for instance [43, 
Theorem 5.1] where a counterexample to uniform exponential stability (or equivalently, to the uni¬ 
form observability estimate (44) without viscosity) is given for a linear damped wave equation. As 
for space semi-discretizations, this is caused by spurious highfrequency modes that appear when 
discretizing in time, and that propagate with a vanishing velocity (as At —> 0). The role of the vis¬ 
cosity term is then to damp out these highfrequency spurious components. Note that other remedies 
to the lack of uniformity are proposed as well in [19, 43], for instance filtering the highfrequencies. 

As before, we focus here on the use of viscosity terms, more appropriate in our nonlinear context. 

The main assumption in Theorem 3 is the uniform observability inequality (44). 

Certainly, the most general result, which can be directly used and adapted in our study, can be 
found in the remarkable article [19] (see also references therein, of which that paper is a far-reaching 
achievement), from which we infer the following typical example of a viscosity operator: 

y At = —(At) 2 A 2 = (At) 2 A* A. 

For this choice, it is indeed proved in [19, Lemma 2.4, and (1.17), (2.17) and (2.20)] that, if 
the observability inequality (9) is valid for the continuous model, then the uniform observability 
inequality (44) holds true for the time semi-discrete model (45). We could take as well the viscosity 
term V A t = — (idjc — (At) 2 A 2 ) -1 (At) 2 A 2 , which yields as well (44), and which has the advantage 
of being bounded. 

Note that, in contrast to space semi-discretizations (see Remark 2), here, for time semi-discretizations, 
the above choice of a viscosity systematically works in order to recover uniform properties. 

Apart from the viscosity term, note that the time discretization is an implicit midpoint rule. 

This choice is relevant for at least two reasons. The first is that an explicit time discretization 
would then immediately lead to a violation of the stability CFL condition, and thus the scheme is 
unstable. Actually, in Section 3.3, we are going to consider full discretizations, and then, to be 
always in accordance with the CFL stability condition, it is better to choose an implicit scheme. 
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This choice is also relevant with respect to the conservation of the energy, for the linear conservative 
equation (45), at least, without the viscosity term. 

Some variants of the midpoint rule and of the design of the viscosity term are possible (see [19, 
Section 2.3]). 

3.2.2 Proof of Theorem 3 

We follow the lines of the proof of Theorem 1. 

As in the continuous setting and in the space semi-discrete setting, the global well-posedness 
of the solutions follows from usual a priori arguments (using energy estimates), and uniqueness 
follows from the assumption that F is locally Lipschitz on bounded sets. 


First step. Comparison of the nonlinear equation (41) with the linear damped model. 

We first compare the nonlinear equation (41) with the linear damped equation with viscosity 


J-fc+l _ gk 

At 

^fc +1 _ ^fc +1 


A 


z k + z k+1 


B 


z k + z k+1 


= 0 , 


(46) 


At 


= V A t z 


fc+i 


Lemma 9. For every solution (u k )k^iN of (41), the solution of (46) such that z° = u° satisfies 

N-l 


E 

k=0 


B 1/2 


z k + ~k+1 


N-l 

^E 

k=0 




B x ' 2 


£||(Va*) 1 / V= +1 ||| + 

At-'"- 1 

=■ E iiw‘ +i i& 

fc=0 


fc=0 

u k + u k+1 \ 

2 N-l 

+ 2E 

X k=0 

N-l 

B l/2 (U k +U k+1 

2 ) 

B F \ 2 


N-l 


+ 2^ \\{V At ) 1 / 2 u k+1 \\ 2 x + At E llVAtu fc+1 |ft- 


k=0 


k=0 


Proof. For every /c, setting ^ k = u k — z k and i/j k = u k — z , we have 


hk+l _ „/,fc / ^pk _|_ ^pk-\- 1 ^ p p yk _|_ yk-\- 1 \ „k _i ^ k 


f) k+1 - %f k 


At 

/,fc+l _ 


= -A 


— B [ F 


(47) 


- - 

At 


Denoting £+ = and taking the scalar product in X in the first equation of (47) with 


il> k+ 1 +il> k ■ 


it follows that 


E^k +1 — -E+ + At (p3 


z K + z 


■fc+i 


Z + Z' 


■k+1 




= -A t(^BF 


u +u 


fc +1 


u + u 


fc +1 


X 


A tpBF 


u k +u k+1 


z k + jfc+i 


+ A t(B 
x \ 


z + z 


- fc +1 


u + u 


fc +1 


< A tfBF 


u + u 


fc +1 


z K + r 


fc+i 


+ At (B 
x \ 


yk+l 


U + U‘ 


fc +1 


X 
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Using the second equation of (47), we infer that 


^ fc+1 -V+ 1 =At||(V At ) 1 /V +1 ||i + ^||V A ^ fe+1 ||^. 


Subtracting the latter equation to the previous one, we obtain that 


E^k+i — E^pk /\t 

< A t^BF 


yk~\~ 1 


2 

u K + u K 


z K + z 


' fc +1 




+ At||(v At ) 1 / 2 v> fc+1 lli + ^l|v At ^ +1 ||^ 


z k + jfc+i 


+ At 
x \ 


z k + z k+1 


u k + u k+1 


Summing from k = 0 to k = N — 1, since ijr = 0, we get that 


N -1 


E,i 


AtJ 2 ( B 


fc=o 


JV-l 


Z fe + z k +1 

2 

JV-l 


z fe + ^fc+i 


\2 JV-l 


+ At ^ ||(V At ) 1/ V +1 l& + ^- E W v ^ k+1 \\ 

k=0 k —0 


^ At ( BF 


u k + u k+1 


z k + z k+1 


X—1 


AtJ2( B 


z k + z k+1 


fc=0 ' v / ' A k -0 

Thanks to the Young inequality, and since E^n ^ 0, we infer that 

JV-l / fc , \ 2 JV-l / A A 2 JV-l 

2 


At 


E 

/c —0 


i ? 1 / 2 


z fe + j-fc+i 


u fc + u fc+i 


X 


X 

X-l 


X - 1 //w \2 iV_1 

At £ ||(v At ) 1 / V + 1 lli- + EE E HWAl 




A*E 


fc—o 


fc =0 

b 1/2 f 


k =0 


+ ~fc+l 


X 


X-l 

A tE 

k—0 


B 1 / 2 


u k + ~k +1 


The lemma follows, using that ||(V At ) 1/2 V' fc+1 ||^ ^ 5 1|(V At ) 1/2 ^ fc+1 113c - ll(V At ) 1/ V +1 |lx and 
that||V At V- fc + 1 ||^^i||V At z fe+1 || 2 f -||V At u fc + 1 ||^. ' ' □ 

Second step. Comparison of the linear damped equation (46) with the time semi-discretized 
conservative linear equation with viscosity (45). 

Lemma 10. For every solution {z k ) k ^iN of (46), the solution ((f> k )k£iN of (45) such that fP = z° 
satisfies 


JV-l 

Y 

k—0 


B 1 ' 2 


kk +1 


JV-l 


At 


JV-l 


+ E ii(v At ) 1 / v + 1 iil- + t E ii W +1 ii 


X k =0 

( X— 1 / u i_i \ 2 X—1 A , X—1 \ 

E Bl ' 2 ( z - Af —) + Eikva,) i/! ^ +i ii 2 y+4^Eiiva,»‘ +i i£j. 

fe = 0 ' ' X k -Q fc _Q J 


k -0 


with kr = max(l + (4T 2 + 1) 2 ||B 1 / 2 || 4 , 2). 
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Proof. Setting 9 l = cj) 1 — z l and 9 l = c/) 1 — z l , we have 


Qi+l _ Qi 

At 

e i+l - e i+1 

At 




Qi + Qi +1 


- B 


y* + l 


= 0 , 


( 48 ) 


= v At e i+1 . 


Taking the scalar product in X in the first equation of (48) with i(0 l+1 + 9 X ), we get 

'z i + z i+1 \ 9 i +9 i+1 \ 


Egi+i — Egi — At { B 




X 


Now, using the second equation in (48), we obtain that 


Eg i+1 - E gi+ i = At ||(V At ) 1/2 ^ +1 Hi + ^tH|Va# + 1 || 2 y . 


(At) 2 


(49) 


(50) 


Subtracting (50) to (49), and then summing from i = 0 to i = k — 1, with k ^ N, using that 
9° =0, we obtain 


Ee- + At X] ||(V At ) 1/2 ^ +1 Hx + \WAt9' +1 \\l 


i =0 


2=0 


k-1 


= AtJ2{B 


i =0 


i + z i+1 \ 9 i + 9 i+1 


(51) 




In passing, note that (51) implies that 


N-1 


zv — i a , N— 1 jv — i I / 

E UVAt) 1/2 e k+1 \\ 2 x + — E IIVA t 0 fc+1 ||^ < E ( B ( 

k—0 k—0 k—0 \ ' 

Now, summing (51) from k = 0 to k = N— 1, using that NAt ^ T, we get 


N-l 


Z /C + Jfc+1\ Qk+Qk+l 


(52) 


N-l 


N-l 


E ^+At e (^ - 1 - fc)ii(VA t ) 1/2 0 fe+1 ii! + ^ E - 1 - fc)iivA t 0 fc+1 n 


X 


k—0 


k =0 


X-l 


= At E^- 1 -*) ( 5 


+£ fc+1 \ Qk + Qk+1 


k—0 \ 

from which it follows in particular that 
N-l N-l 


X 


k—0 
N-l / 

k—0 \ 


k _|_ jfc+1 \ Qk _(_ $fc+l 


E^< T E B 


fc =0 


fc =0 


+ Jfe+l \ Qk + flfc+l 


( 53 ) 


X 


Besides, thanks to (50), since 0° = 0, we have 


N—2 


N-l 


N—2 


N—2 


IV —^ IV —JL IV —4 / A |\2 iV —^ 

E ^ + 1 = E E ll(VA t ) 1/2 0 fc+1 ||i + V 1 E l|VA^ fc+1 ||3o 


k—0 


k—0 


k—0 


k—0 
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and, using (52), we infer that 


N—2 


N—l 


N-1 


fc =0 


y! ^/c+i ^ E 0 k +At ^ ^ 5 

/c —0 /c—o 

Using (49) for i = iV — 1, we have 

Eg N = EgN- 1 + At ( B 


+ 2 fc+1 \ 9 k + 8 k+1 \ 


/ 5 

/ X 


v iV- 


_1 + 2 ' JV \ 


7 X 


(54) 


(55) 


Summing (54) and (55), we obtain 


N-1 


N—l 


N—l 


y -E'0fc+i ^ 2 ^ Egk + 2At ^ 5 


+ z k+1 \ e k + 6> fc+1 \ 


fc =0 


/c=0 


k—0 


! x 


Finally, summing three times (53) and (56), we get 


(56) 


N—l 


N-l 


£(£ 9 *+£fc> + 1 )<4T£ (B 


+ 2 fc+1 \ e k + 6> fc+1 \ 


Noting that 


/c— 0 

0 fc + ^fc+i 

2 


/c— 0 


/ X 


(57) 


< 5 (ll^llx + ll^' +1 Hx) = E gk +Eg k+1 , we infer from (52) and (56) that 


N-l 

E 

k—0 


e k + e k+1 


N-l 


At 


N-l 


+ £ ||(v A t) 1/2 ^ +1 lli- + ^ £ IIVA t 0 fc+1 ll 


X 


fc =0 

N-l 


i(iT + l)J2{B' + et + dk+ 7 


k—0 


1 . 


(58) 


By the Young inequality, we have 


AT-l 


(4T+1)£ (b 


k =0 


7 , 


X 
N-l 

<-(4T+l) 2 || J B 1 / 2 || 2 £ 

fc =0 


B 1 / 2 


+ z k+1 


2 JV-1 

+^E 

* k—0 


Qk + Qk +1 


and therefore we get from (58) that 


N-l 

-Y 
2 ^ 

k=0 


+ 6> fc+1 


N-l 


At 


N-l 


Y ll(VA t ) 1/2 ^ +1 |lA + T E \\VAtO k+1 \\ 


X k—0 


k—0 


<i(4T+l) 2 ||B 1 /2||2^ 
k—0 


N-l 


B 1 ' 2 


yk+1 


(59) 


X 
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Now, since <jr = 0 k + z k and (\r = Q k + z k , using the inequality (a + b ) 2 ^ 2(a 2 + fr 2 ), we have 


N—l 

E 

fc =0 


B 1 / 2 


ife , jfc+i 


2 

JV-l 


2 JV-l 


At 


JV-l 


+ E ||(v At ) 1 / V +1 lll + ^ E HW +1 ILl 


x fc =0 


/c —0 


< Iis 1/ T £ 

k=0 
N—l 

+ E 


+ # fc+1 


AT—1 


N-l 


II- 


/c —0 


B 1//2 f 


z fe + z k+1 




+ 2^ ||(V At ) 1/2 0 fc+1 ||^+At £ ||V At 0 fe+1 ||) 

x fe =0 fc=o 

2 JV-l JV-l 

+ 2£ ||(V At )VV* 1 |& + Ai £ ||V At . fe+1 |||, 




fe=o 


/c=0 


and, using (59), the lemma follows. 


□ 


Third step. Nonlinear energy estimate. 

We define w by (17), as before, with /3 > 0 to be chosen later. 

Lemma 11. For every solution (u k )k^iN of (41), we have 


JV-l 


At w(E c / > o) 


k=0 


< 


B 1/2 


U k + M fe+1 


B l /2 F ( u k +U k +' 


X , 


JV-l 


||J3||ATAt H*(w(E ( / ) o)) + (w(F» + l)At ^ (# 


fc =0 


M fc + U k+1 


,F 


u k + u k+1 


Proof. For every k € IN, we set f k = U 1 u k and f k = U 1 u k . By definition of p, we have 
p(/ fc ) = U~ l F{Uf k ) for every k £ IN. We set £o = min(l, p(so)) and we define, for every k £ IN, 

the set flj = {x £ ft | | - + ;f + (x)| < £o}- Using (6), we have p (^—- 5 —— j 6 [0,s 2 ] and 


1 //* + / 


fc 1 ffc+l 


fa*, b dp jQk 4 


bdp G [0, s 2 ]. 


Since 17 is convex on [0,Sg], applying the Jensen inequality, we get 


H 


l I f k _|_ rfe+i 


./nf cl 


T 2 P 


bdp 


€ 


f k + / fc+1 


f Q k bdp J n k c 2 
Using (6) and the sign condition (5), we infer that 


<?(- 

\ C 2 


f k + f k+1 


b dp. 


H 


1 




fk _|_ yfc+i / fk _|_ y^+l 


JnkbdpJn* 4 

Since U is increasing, we deduce that 
2 


C 2 f^k bdp J Qfc 2 


dp. 


f k + f k+1 


inf 


b dp ^ 4 I bdp Ft 1 [ — -jr —[ b- 

jQk y c 2 ^ d[J, •/fif 


fk _|_ y/c+1 / fk _|_ yfc+1 \ 


2 J d "J' 
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and therefore. 


2 


N-l . 

At y w(E (/) o) / p 
k=o J n f 


fk _|_ yfc+i 


bdp 


n i t- (lit 

< At V Cn / bdp w{E 1 ko)H~ 1 - T ——— / b 

tr'o Aif \ c 2 Jnibdp J n $ 


fk _|_ y/c+l / fk _|_ y/c+1 



dp . 


Hence, according to the Young inequality ^ H(A) + H*{B) : we get that 


p/c _l_ f/c+l / j?k _|_ yfc+1 


y /* / £k _i_ f/c+i \ - y 7 y /* f^-L f^+i 

Al^w(V) / P-o- &d/i^c 2 Aty / b- --- p 

k—0 dn i V / fe=0 

+ NAtc% [ bdp H*(w(E r )). (60) 

Jn 

, / jrfc_i_77e + l\ ffc | 77c+l 

Besides, in £7 \ fif, using ( 6 ) we have p ( 7 ^ < 7 ^ - ? and using (5), it follows that 


'f k + f k+1 ' 


/ n\nj 


2 


fk _|_ y/c+1 / fk _|_ y*&+l 


dp. 


Using this inequality and (60), we obtain 


N-l 


At y w(E c/> o) / p 


fk _|_ fk+1 


k—0 


bdp<NAt / bdp H*(w(E ( f > o)) 
J n 


iv-i 


J ’_y /• fk I ck+l / fk I rfc+1 \ 

+ w*>) + i)At y / b J y P 7 / d/x. 


fc =0 


7fc + l 

Now, defining as in Lemma 3, we set = {x £ fl | | J + .{ -(x)| ^ £i}. Using ( 6 ), we have 


r+r 


S [ 0 , s§] and 


1 


/nj bdp J^k 


/' 

fk + fk+1 


2 


bdp £ [ 0 , Sq]- 


Since H is convex on [0,Sg], by the Jensen inequality, we get 
H 


( i f 

f k + fk+i 

2 ) 

. 1 [ 

f k + Jk+i 

s ( 

fk + fk+i 

\[ n JjdpJ n k 

2 


Jnj b d d 7n§ 

2 

2 


€ 


1 1 


fk _|_ / fk _|_ y*&+i 


Cl U bdp J n 


bdp 


dp. 


Since H is increasing on [0, Sq], we deduce that 


iv-i ,, 

At y w(E c/) o) / 

fc =0 < 


/' 

f k + jk+i 


2 


b dp 


< At y w^Effro ) [ bdp H 1 ( — r 1 [ b 

y, Jn 2 fc \ c iJn^ bd ^Jn 


fk _|_ y/c+1 / fk _|_ yfc+1 



dp . 
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It then follows from the Young inequality AB ^ H(A) + H*(B) that 


N -1 


At w(E t / > o) j ^ / — bdp ^ NAt f bdp, H*(w(E t p o)) 

k^o 2 -' n 


N -1 


Cl 


At E 


6 / fc + / fc+1 p / / fc + / fc+1 


k =0 


The estimate in £l\£l k is obtained similarly. The lemma is proved. 

Fourth step. End of the proof. 

Lemma 12. We have 

E U N < E u 0 ^1 - ptL~ x (^fj- 

for some sufficiently small positive constant pt ■ 

Proof. Summing in k the energy dissipation relations (43), we get 


dp. 


□ 


(61) 


N—l 


E u n - E u o = -At ( B 


k =o 


u k + u k+1 


,F 


u k + u k+1 


N-l 


-A tJ2 ll(VA t ) 1/2 ^ +1 Hx - ^r- E \\VAtu k+1 \\ 2 x- (62) 


2 N-l 


k -0 


k =0 


Using successively the observability inequality (44) and the estimates obtained in Lemmas 10, 9, 
and 11 , and since NAt ^ T, we get 


Ctw(E^o)U°\\ 2 x 


N-l 


< 


k T \\B\\TH*{w{E r )) + k T {w{E r ) + l)Ai E ( B 


k =0 

/N-l 


U k + u k + l 




u k + U kJrl 


N-l 


+ kr&twiEjo) ^ \\(V&t) l/2 u k+1 \\ 2 x + At Y2 \\V^tu k+1 \\ 


, k —0 


k —0 


Recalling that E^o = E u o, using (62), we infer that 

C t w{E u o)E u o < k T \\B\\TH*(w{E u o)) + k T {w{E u o) + 1) (E u o - E u n) . 
Then, reasoning as in the proof of Lemma 4, we get the conclusion (we skip the details). 


□ 


Let now p £ IN be arbitrary and k £ {0,... ,N}. The sequence (v k )k^{o,...,N-i} defined by 
v k = u k+pN for k £ {0,..., TV — 1}, satisfies (41) and thus (61). Noting then that E v n = E u ( P +i)n 
and E v o = E uP n , we deduce that Ek+i — Ek + prM(Ek) ^ 0, where M(x) = xL~ 1 {x) for every 
x £ [0, s§] and E k = E u kN IP- 

As at the end of Section 2.3, setting K r {r) = f* we have 


M(E p ) ^ 


1 

— mm 
Pt f€{o,...,p} 


K r (pr(p-l)) 

£+1 
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We set t = pT. For any 9 £ (0,t], we set l = [|r] £ {0,... ,p}. We have 


E p < M _1 


— inf 

PT 0 <6^.t 




As in the proof of [5, Theorem 2.1], we deduce that 


E uP n A j3L 


1 


Ip 1 ( p T P) 


(63) 


for every p ^ 1/(ptH'(sq )), with ip defined by ( 8 ). 

Moreover, under (11), we get that E uP n A j3 (7 / 7 ) -1 f° r every p sufficiently large, for a 

certain 73 > 0 not depending on p, At, Eq. Let now k > T / At be a given integer. We set p = k/N. 
Since pN A k and thanks to the dissipation property (43), we have E u k A E uP n. Besides, we have 
PtP > Pr(k/N — 1). Since ip and L are nondecreasing and thanks to (63), it follows that 


E uk A QL ^ ( kAf _ , 

for every k ^ pt h’(s 2 ) ~tt' Moreover, under (11), we have E u k A /3(7F ) _1 for k 

sufhciently large. Theorem 3 is proved. 


3.3 Full discretization 

Following [19] (see also references therein), results for full discretization schemes may be obtained 
from the previous time discretization and space discretization results, as follows: it suffices to 
notice that the results for time semi-discrete approximation schemes are actually valid for a class 
of abstract systems depending on a parameter, uniformly with respect to this parameter that 
is typically the space mesh parameter Ax. Then, using the results obtained for space semi- 
discretizations, we infer the desired uniform properties for fully discrete schemes. 

More precisely, the class of abstract systems that we consider is defined as follows. Let ho, B, T\, 
72 , Ci, C 2 , 1C, 17, za> and v 3 be positive real numbers, let so £ (0, 1 ], and let g : 1R —> R be a function. 
We define ^(h 0 , B, 71, 72 , Ci, C 2 , /C, g, s 0 , 17 , i/ 2 , ^ 3 ) as the set of 5-tuples (,i?7, ,e4- D{.sA h ),AB h , &h), 
where, for every h £ [0, ho): 

• SCh is a Hilbert space (of finite or infinite dimension), endowed with the norm || ||^; 

• s^h '■ C 3&h —t is a densely defined skew-adjoint operator; 

• ABh : SAh —t SAh is a bounded selfadjoint nonnegative operator such that \\ABh\\ A E; 

• there exist £A h £ [ 71 , 72 ] and % G [Ci,C 2 ] such that 

VhWhml < J* h (\\K /2 Mt)\\i + h°\\rf /2 Mt)\\i) dt, 

for every solution of 

+^h(ph{t) = 0 ; 

• : 3Ah —> AAh is a mapping that is Lipschitz continuous on bounded subsets of 3Ah> with 
Lipschitz constant less than 1C, and satisfying (5) and ( 6 ) with the function g\ 
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• g is an increasing odd function of class C 1 such that g( 0) = </(0) = 0, lim s _,.o sg'(s) 2 /g(s) = 0, 
and such that s K > -\/sg(\/s) is strictly convex on [ 0 , Sq]; 

• Yh : 5£h —> 9£h is a positive selfadjoint operator, and the family (ft CT / 2 \\V^ 2] \)he{o,ho) 
uniformly bounded; 


• any solution Uh of 

u' h {t) + £/ h u h {t) + @h&h{u h {t)) + h a y h u h {t ) = 0, 


with Uh{ 0) £ D{e/ h ), is well defined on [0, +oo), and satisfies 
\\uh{t)\\l < y 3 max(^, \\u h {0)\\ h )L 

for every time t ^ 0, where L and if are defined (in function of g) by (7) and ( 8 ). 

Within the framework and notations introduced in Sections 1 and 3, under the assumptions of 
Theorems 1 and 2, there exist positive real numbers B, 71, 72 , Ci, C 2 , /C, V 2 and 1 / 3 , such that the 
5-tuple (X, A, D(A), B, F ) and the one-parameter family of 5-tuples A& x , Xa x , Ba x , F/\ x ), 

A x £ (0, Ax 0 ), belong to the class ^(Axo, B, 71, 72,Ci,C 2 , /C, g, s 0 , vi, zz 2 , ^ 3 )- 

Here, the parameter ft used in the definition of the abstract class above stands for the space 
semi-discretization parameter Ax, whenever ft > 0, and if ft = 0 then we recover exactly the 
continuous setting of Section 1. 

We claim that Theorem 3 can be applied within this class, uniformly with respect to the 
parameter ft. This can be checked straightforwardly, noticing in particular that the constants 
appearing in the estimates of that result depend only on the data defining the class. 

In particular, from that remark, we infer full discretization schemes of (1), in which we have 
first discretized in space, and then in time, obtaining an energy decay as in Theorems 1, 2 and 
3, uniformly with respect to the discretization parameters Ax and At. Said in other words, we 
apply Theorem 3 to the one-parameter family of systems given by (28), parametrized by Ax, and 
for which, by Theorem 2, we already know that the solutions decay in a uniform way. We have 
thus obtained the following result. 


Theorem 4. Under the assumptions of Theorems 1, 2 and 3, the solutions of 


^ 7 fc + 1 _ nj‘ 
u Ax u 


Ax 


At 


U Ax U Ax 


■ A/\ x 


‘Ai 


yk+1 ' 
l Ax 


+ BaxF/\x 


,k 1 yk+1 ' 

l Ax + a Ax 


VAx 


‘Ai 


y/c +1 ' 
l Ax 


= 0 , 


= Va t u k + x \ 


At 


are well defined for every integer k, for every initial condition u° Ax £ X/± x , and for every Ax £ 
(0, Axq), and the energy of any solution satisfies 


< Tmax( 7 l ,£ BL )ft ^_ 1( ^ Af ) j , 

for every integer k, with 72 ~ Ct/T{ l+e 2T ^ B ^ max(l, T||H||)) and 71 ~ ||H||/ 72 . Moreover, under 
(11), we have the simplified decay rate 


E„ 


< Tmax^iy^) (#') 1 



for every integer k, for some positive constant 73 ~ 1. 
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Example 1. Let us consider the nonlinear damped wave equation studied in Section 2.2.2. We 
make all assumptions mentioned in that section. 

We first semi-discretize it in space by means of finite differences, as in Remark 2, with the 
viscosity operator V Ax = — (Ax) 2 Aak- At some point x a of the mesh (<r being an index for the 
mesh), we denote by u AXta (t) the point of IR” representing the approximation of u(t,x a ). Then 
the space semi-discrete scheme is given by 


u Ax,cr(t) A Ax U/\. Xt< 7 {t) b Ax {x a )p{x a ,U Ax a (t)) (Aa;) — 0 ) 

and the solutions of that system have the uniform energy decay rate L(l/'i/> _1 (t)) (up to some 
constants). 

Now, discretizing in time, we obtain the numerical scheme 


_ 7/ fc _i_ 7 j k + 1 

^ Ax,cr “'Ax,cr “Ax,a ' “Ax,<r 

At ~ 2 : 

~fc+i _ k 
^A x,a ^A x,a 


At 


+7, k+l ( v k +v k+1 ' 

A ^ Ax ,<j ^ ^ Ax,a .7 / \ / “Ax,a ' “Ax,a 

x,o 7) “Ax,cr\pC<j)P I *^<7■ 


~fc+l 


- (Ax) 2 A AXt0 


~k+i _ fc+i 
U Ax,a a Ax,a 


- = -(A i) a AL,^ 


At 


~fc+i _..fe+i 

^Ax,a ^ A x,o 

At 


= -{AtfAl^v^ 


= 0, 


and according to Theorem 4, the solutions of that system have the uniform energy decay rate 
L(l/'i/> _1 (t)) (up to some constants). 


4 Conclusion and perspectives 

We have established sharp energy decay results for a large class of first-order nonlinear damped 
systems, and we have then studied semi-discretized versions of such systems, first separately in 
space and in time, and then as a consequence, for full discretizations. Our results state a uniform 
energy decay property for the solutions, the uniformity being with respect to the discretization 
parameters. This uniform property is obtained thanks to the introduction in the numerical schemes 
of appropriate viscosity terms. 

Our results are very general and cover a wide range of possible applications, as overviewed in 
Section 2.2. 

Now several questions are open, that we list and comment hereafterin. 

(Un)Boundedness of B. The operator B in (1) has been assumed to be bounded, and this 
assumption has been used repeatedly in our proofs. This involves the case of local or nonlocal 
internal dampings, but this does not cover, for instance, the case of boundary dampings. To give 
an example, let us consider the linear ID wave equation with boundary damping 

f d tt u-d xx u = 0, t € (0, +oo), x € (0,1), 

\ u(t, 0) = 0, d x u(t, 1) = -ad t u(t, 1), f € (0, +oo), 

for some a > 0. It is well known that the energy E(t) = | /' (\d t u(t, x)\ 2 + \d x u(t, x)| 2 ) dx of any 
solution decays exponentially (see, e.g., [12]). It is proved in [41] that the solutions of the regular 
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finite-difference space semi-discrete model with viscosity 

f ^ (Ax) ^Aa:^i x — 0, 

^Ax(f; 0) — 0; D^ x u/\ x (t : 1) — ad t u(t, 1), 

where D& x is the usual ID forward finite-difference operator, given by 

-1 1 ... 0 \ 

0 : 

: 1 

0 ... 0 - 1 / 


° Ax ~ Ax 


have a uniform exponential energy decay. Without the viscosity term, the decay is not uniform. 
This example is not covered by our results, since the operator B in that example is unbounded. 

We mention also the earlier work [11], in which the lack of uniformity had been numerically 
put in evidence for general space semi-discretizations of multi-D wave equations with boundary 
damping. By the way, the authors of that paper proposed a quite technical sufficient condition 
(based on energy considerations and not on viscosities) ensuring uniformity, and applied it to mixed 
finite elements in ID and to polynomial Galerkin approximations in hypercubes. It is not clear if 
such considerations may be extended to our nonlinear setting. 


More general nonlinear models. In relationship with the previous problem on B unbounded, 
the question is open to treat equations like (1), but where the nonlinearity F involves as well an 
unbounded operator, like for instance the equation 

u'(t) + Au(t) + BF(u(t),Vu(t)) = 0. 

There, the situation seems widely open. In our approach, what is particularly unclear is how to 
extend Lemma 3. 

An intermediate class of problems that has not been investigated at the discrete level is for 
instance the class of semilinear wave equations with strong damping 

dttu — Au — aAdtu + bdtu + g * A u + f{u) = 0, 

with / being not too much superlinear. Here also, many variants are possible, with boundary 
damping (see [5, 8], with nonlocal terms (such as convolution), etc. There exists a huge number 
of papers establishing decay rate results for such equations, see e.g. [15, 20] and the references 
therein and see also [7] for a nontrivial extension of the optimal-weight convexity method to the 
case of memory dissipation (non-local dissipation) but nothing has been done for discretizations. 

Geometric conditions and microlocal issues. Another class of equations of interest, not 
covered by our main result, is the stability of semilinear wave equations with locally distributed 
damping 

d u u — Au + a(x)d t u + f(u) = 0, 

with Dirichlet boundary conditions. Here, a is a nonnegative bounded function assumed to be 
positive on an open subset oj of f2, and the function / is of class C 1 , satisfying /(0) = 0, s/(s) ^ 0 
for every s £ 1R (defocusing case), |/'(s)| < 1 with p < n/(n — 2) (energy subcritical). We 

set F(s) = fg f. It is proved in [ | (see also some extensions in [16, 22] and a variant in [10]) that, 

under geometric conditions on w, the energy (which involves here an additional nonlinear term) 

J n (\( dtU ^ + ^ll Vw H 2 + F ( U )) dx 
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decays exponentially in time along any solution. It is natural to expect that this exponential decay 
is kept in a uniform way for discrete models, if one adds appropriate viscosity terms as we have 
done in this paper. Note that this is not covered by our results since the above energy involves an 
additional nonlinear term. 

Besides, in the general case where lo satisfies the Geometric Control Condition (GCC) of [12], 
the arguments of [10, 16, 22] rely on microlocal issues and it is not clear whether or not such 
arguments may withstand discretizations. For instance, it is not clear what GCC becomes in a 
discrete setting. It is also not clear what a microlocal argument is at the discrete level, and such 
considerations may lead to several possible interpretations. In brief, we raise here the completely 
open (and deliberately informal and imprecise) question: 

Do microlocalization and discretization commute? 

Uniform polynomial energy decay for linear equations without observability property 

Let us focus on linear equations, that is, let us assume that F = idjf in (1). It is well known 
that, in the continuous setting, the observability inequality (9) holds true for all solutions of the 
conservative linear equation (15), if and only if the solutions of the linear damped equation (13) 
(which coincides with (1) in that case) have an exponential decay (see [21]). 

If the observability inequality (9) is not satisfied, then the decay of the energy cannot be 
exponential, however, it may be polynomial in some cases. It may be so for instance for some 
weakly damped wave equations in the absence of geometric control condition (see [12]) but also for 
indirect stabilization for coupled systems, that is when certain equations are not directly stabilized, 
even though the usual geometric conditions are satisfied (see [2, 3, 4]). In that case, it would be 
of interest to establish a uniform polynomial decay rate for space and/or time semi-discrete and 
full discrete approximations of (1). In [1], such results are stated for second-order linear equations 
(certain examples being taken from [3, 4]) , with appropriate viscosity terms, and under adequate 
spectral gap conditions. 

Extending this kind of result to a more general framework (weaker assumptions, full discretiza¬ 
tions), and to our nonlinear setting, is open. Our strategy of proof is indeed strongly based on the 
use of observability inequalities. 
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